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ORTHOGONAL AND SYMPLECTIC n-LEVEL DENSITIES 


A.M. MASON AND N.C. SNAITH 


Abstract. In this paper we apply to the zeros of families of L-functions with orthogonal or sym- 
plectic symmetry the method that Conrey and Snaith m used to calculate the n-correlation of 
the zeros of the Riemann zeta function. This method uses the Ratios Conjectures m for averages 
of ratios of zeta or L-functions. Katz and Sarnak m conjecture that the zero statistics of families 
of L-functions have an underlying symmetry relating to one of the classical compact groups U{N), 
0{N) and USp{2N). Here we complete the work already done with U{N) to show how new 
methods for calculating the n-level densities of eigenangles of random orthogonal or symplectic ma¬ 
trices can be used to create explicit conjectures for the n-level densities of zeros of L-functions with 
orthogonal or symplectic symmetry, including all the lower order terms. We show how the method 
used here results in formulae that are easily modified when the test function used has a restricted 
range of support, and this will facilitate comparison with rigorous number theoretic n-level density 
results. 
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1. Introduction 

1.1. Motivation. Long-standing results in random matrix theory show that the n-point correlation 
functions (in some situations these are also called n-level densities) of eigenvalues from ensembles 
of random matrices such as U{N), SO{N) and USp{2N) can be written concisely as n-dimensional 
determinants of matrices whose elements are the kernel belonging to the particular ensemble (see, 
for example, [16] where the kernels and correlation functions are written down for these groups). 
These results are elegant and exact, and their determinantal form is very useful for calculations 
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within random matrix theory (RMT). However, the corresponding quantity in number theory, the 
n-point correlation function (or level density) of the complex zeros of an L-function, or of families of 
L-functions, does not seem to take this concise determinantal form, except in a suitable asymptotic 
limit where the statistics are expected to agree with RMT (see da [Ml EH] for reviews on the 
connection between number theory and random matrix theory). This leads us to a question: is there 
a different, albeit less elegant, way to write the random matrix eigenvalue correlation functions that 
might help to make precise conjectures about the form of the correlations of zeros? Furthermore, 
given that the determinantal form is not available in the number theory case, is there a natural 
form in which to write the n-correlations of the zeros that is useful for further applications? 

The first question has been answered in the case of the Riemann zeta function, which has 
zeros displaying the statistics of eigenvalues of matrices from U{N), with Haar measure, in the 
appropriate limit. The n-point correlation functions of the Riemann zeros, including the lower order 
correction terms to the limiting random matrix theory result, were first studied by Bogomolny and 
Keating PE] and then by Conrey and Snaith [291 EZj- In particular, Conrey and Snaith derive 
the n-point correlation function of U{N) eigenvalues using average values of ratios of characteristic 
polynomials in a method that is precisely reproducible in the number theory case using a conjecture 
on average values of ratios of the Riemann zeta function (see [21] and also [26]). 

It is this first question in the case of matrices from SO{2N) and from USp{2N) and the families of 
L-functions associated with them that is the main focus of this current paper. The n-level density 
functions of eigenvalues from matrices in SO{2N) are calculated in Section [2] using averages of 
ratios of characteristic polynomials, with the final result appearing in Theorem 12.111 The steps 
are mirrored in Section |T| using the Ratios Conjecture for a family of L-functions associated with 
elliptic curves. The zeros of this family are expected to behave statistically like the eigenvalues of 
SO{2N). The Ratios Conjecture is used to derive the n-level desities of zeros of the L-functions 
in this family, complete with lower order terms. This result can be found in Theorem 15.111 In 
Section [3| eigenvalues of matrices from USp{2N) are considered in the same way, see Theorem 13.71 
and Section [6| demonstrates the method with a family of L-functions showing symplectic symmetry. 
The result can be found in Theorem 16.61 

As reflected in the second question, the motivation for this work is not simply to derive the 
detailed conjectures for the n-correlation of zeros of families of L-functions, but to obtain formulae 
that are useful for further applications. For the past 40 years, since Montgomery’s pioneering work 
on the statistics of the two-point correlation function of the Riemann zeros m, there has been a 
catalogue of work where a rigorous result on n-correlation functions (often called n-level densities in 
number theory literature) of the zeros of an L-function or family of L-functions is derived and then 
an attempt is made to match it with either the limiting RMT correlation function (for example 
[sniiiaiiMiiHniiMiEi]), or a conjecture derived using a Ratios Conjecture that includes lower order 
terms (for example [Ml EHl HB [TOl Il9]l. For general n this comparison is difficult, and is made 
more so by the fact that all the rigorous number theoretical results hold only for test functions with 
restricted support (of their Fourier transform) and by the fact that if comparison is made with the 
determinantal RMT form, this has nothing like the shape of the number theoretical results. By 
test function we mean the function / in the definition of the n-level density ()5.ip , which is sampled 
at the positions of an n-tuple of zeros. There is an equivalent test function in the definition of the 
n-level density of eigenvalues, for example in (12.11) . that is sampled at the positions of an n-tuple 
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of eigenvalues. In [28] Conrey and Snaith demonstrate that in the U{N) random matrix theory 
case, the form of the n-point correlation function resulting from the method using the average of 
ratios of characteristic polynomials allows for immediate simplification when the support of the test 
function is restricted - something that is certainly not true of the determinantal form. This allows 
them to apply to the re-point correlation function for eigenvalues from U(N) the same restriction to 
the support of the test function that was used by Rudnick and Sarnak in [80] when looking at the 
re-point correlation function for zeros of a general L-function. Then the identical structure of the 
two expressions reveals that they coincide in their respective asymptotic limits. This addresses the 
second question above in the case of the zeros of the Riemann zeta function. The second question 
in the case of SO{2N) or USp{2N) matrices is approached in Section [7| where it is shown that 
the form of the random matrix theory re-level density given in this paper leads easily to a simpler 
result when the support of the test function is restricted. It will be the subject of future work to 
show whether this allows more straightforward confirmation of the limiting random matrix form 
for rigorous calculations of the re-level density for associated families of L-functions. 


1.2. Background to the Ratios Conjectures. The connection between RMT and number the¬ 
ory became apparent in the 1970s when Montgomery after a conversation with Dyson [32], 
conjectured that in the appropriate scaling limit the pair correlation of the zeros of the Riemann 
zeta function have the same form as the equivalent statistic for eigenvalues of unitary matrices 
of large dimension drawn at random from the group U{N) endowed with Haar measure. There 
followed work, both of an analytical and a numerical nature, that provided more evidence that 
zeros of the Riemann zeta function high on the critical line and eigenvalues of large random unitary 
matrices show the same local statistics when both sets of points are scaled to have a mean spacing 
of unity (see for example [SEiiisiiTSiiin]). Even decades later, Montgomery’s conjecture is proven 
only for test functions with restricted support. 

Katz and Sarnak [561 [57] then considered statistics of zeros of L-functions near the critical point 
(the point at which the real axis crosses the critical line on which the complex zeros are expected to 
lie). They predict that in a natural family of L-functions these zeros would, when averaged across 
the family, display the same statistical behaviour as the eigenvalues near 1 of matrices chosen at 
random with respect to Haar measure from the classical compact groups U{N), 0{N) or USp{2N). 
As in the case of the Riemann zeta function, this correspondence is expected to occur when the zeros 
and eigenvalues being considered are scaled to have a mean spacing of unity and in an appropriate 
asymptotic limit: normally the zero statistics tend to the eigenvalue statistics of large matrices as 
a parameter intrinsic to the definition of the family tends to infinity. The correspondence with 
random matrix theory has been proven for specific families and for test functions with restricted 
support (see for example [111 Ell i39i [ill iMi ia [53i Eg [nl Ei 1761 [771 iig [89^ 

It is a long-standing conjecture that the moments of the Riemann zeta function grow asymptot¬ 
ically like 


1 

T 



|C(l/2 + R)p^dt 


rsj 
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for large T, where we have chosen to write the coefficient in this way for ease of comparison later 
on with random matrix theory. Here ax is a product over primes 



but gx is very difficult to calculate. Moments of the Riemann zeta function have been the subject 
of study for one hundred years, but still gi = 1 [l3] and <72 = 2 m are the only proven values. 
Conjectures have risen from number theoretical methods, for example [24112511^ . but at the start 
of this millennium there was no viable conjecture for any values beyond k = 4 and random matrix 
theory provided the only available method for predicting higher moments. 

The characteristic polynomial, Ayi(s), of a random unitary matrix. A, is used to model the 
Riemann zeta function. The analogous quantity to (II.ip is 


(1-3) Mn{X) = f |AA(l)P'*'(iAHaar, 

JU(N) 

where the integration is with respect to Haar measure. This average can be computed using 
Selberg’s integral (see [66] ), giving 


(1.4) 


Mat (A) 


^ r{j)r{j + 2X) 
(r(j + A))2 


where the final line is the asymptotic for large N. The Barnes double gamma function |1| can be 
used to express /(A) concisely: 


(1.5) 


/(A) 


(g(l + A))^ 

G(1 + 2A) 


It is conjectured [60] that 

( 1 . 6 ) 


/(A) 


gx 

r(i + A2)- 


This conjecture agrees with established or independently conjectured results for gk for k 


1 , 2 , 3 , 4 . 


We note that when we compare (|1.1D and (|1.4p we see N is playing the role of log T, an equivalence 
that is also observed by equating the density of zeta zeros to the density of eigenvalues. This 
observation is key, allowing comparison of random matrix results to number theoretic quantities 
even when the limits N ^ 00 and T —)■ 00 have not been reached. 


For a family of L-functions, a similar quantity is the average, over the family, of the values of 
the L-functions at the critical point. If we denote by L” a family of L-functions Lf{s), indexed by 
/, and is the number of members of the family, then the moment 


E(%(V2))" 


(1.7) 
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can be modelled by a quantity similar to (jl.3p . only with the average performed over the appropriate 
subgroup of U{N) comprised of unitary, orthogonal or symplectic matrices [I7l|59]: e.g. 

(1.8) / {AAil))^dAHaar, A G SO(2N). 

JS0(2N) 

It was described in |19j that more precise conjectures for the moments of zeta and L-functions, 
including lower-order terms, can be obtained from the “shifted” moments. These are quantities of 
the form 


C(l/2 ~\~ it -\- cti) • • • ^(1/2 -\- it -\- OLk) 

(1.9) xC(l/2 -it- ak+i) ■ ■ ■ C(l/2 -it- a 2 k)dt. 

In that paper Conrey, Farmer, Keating, Rubinstein and Snaith propose a “recipe” for conjecturing 
a precise expression for such moments and for similar averages over families of L-functions. One 
argument in support of the validity of these conjectures, is that the expression produced has an 
identical structure the analogous (and rigourous) moment calculation in random matrix theory for 
a quantity such as m (see also P): 



(1.10) / AA(e-“i)...AA(e-“'=)AA*(e“''+i)AA.(e“2'=)d^Haar. 

JU(N) 

Here A* is the conjugate transpose of the matrix A G U{N). 

A further generalization was made by Conrey, Farmer and Zirnbauer to averages of ratios of the 
Riemann zeta function, 

nf=i C(l/2 + it + ak) Ui=K+i C(l/2 -it- at) 


( 1 . 11 ) 


1 

f 


n?=i C(l/2 + it + 7q) nil C(l/2 -it- 5 r) 
and to averages of ratios of L-functions 

1 ^ nti Lf {1/2+ ak) 


-dt. 


( 1 . 12 ) 




E 


f£jr nil ^f{l/2 + 7q) 

In |21j they develop conjectures for these ratios that are informed by analogous random matrix 
ratio calculations [SollIllES]. 

The “recipe” for creating shifted moment or ratio conjectures involves several steps where we 
neglect terms that may be of the same size as the main term that emerges at the end of the 
procedure. The miracle is that these neglected terms appear to cancel out. The recipe appears to 
produce a correct conjecture down the the level of the constant term or below. 

The main steps of the recipe are as follows; these are reproduced from m, although somewhat 
simplified as we don’t need their full generality for the families with orthogonal or symplectic 
symmetry considered here. For particular cases where a Ratios Conjecture is worked though, see 
for example [26] or |48j . 
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The recipe: Here we assume we have a family T of L-functions, indexed by / E Each 
L-function in the family has a functional equation of the form (assuming we have normalised so 
that the critical line has real part 1/2) 

(1.13) £/(s)=e/T;(s)Z7(l-s). 

As described in [19], we measure the “size” of an L-function by c(/) = |A’j(^)|. This is the logarithm 
of the quantity which for certain families of L-functions is commonly referred to as the “conductor”, 
and by the argument principle c(/)/(27r) gives the density of zeroes near the critical point. 

The L-function also has an approximate functional equation of the form 


(1.14) Lf[s) = ^ + e/T/(s) + remainder, 

for appropriate ranges of summation over n. In addition, for each L-function we can write its 
reciprocal as a Dirichlet series 


(1.15) 


1 ^ ^ L/(^) 


for appropriate coefficients /i/(n); for example, in the case of the Riemann zeta function the Dirichlet 
series for the reciprocal of zeta features /r(u), the Mobius function, which is multiplicative and is 
equal to -1 when n = p is prime and is equal to 0 when n = with i> 


We are interested in understanding the ratio 


(1.16) 


Cf{s]OLK]lQ) 


Cf{s -I- ai) • • • Cf{s + Ok) 
^f{s+ -/!)■■■ Cfis + 'yq) ’ 


averaged over / € L”. For the L-functions we encounter in this paper, L/(s) = Cf{s). 

The first step of the recipe is to replace each L-function in the numerator with the two terms from 
its approximate functional equation ()1.14l) . leaving out the remainder term. Then each L-function 
in the denominator is replaced by the series (|1.15p . Multiplying out the K approximate functional 
equations in the numerator results in 2^ terms. Each of those terms can be written in the form 


(1.17) 


(product of ej-factors)(product of Ay factors) (summand). 

ni,...,nK+Q 


In this expression the recipe indicates we replace each product of e/’s by its expected value when 
averaged over the family. Similarly, the summand is replaced by its average over the family. The 
summations that resulted from the approximate functional equation are over a restricted range of 
integers. The next step is to extend the sums to the full range from 1 to infinity. If the result of 
these steps is denoted as M/(s; Tg), then the corresponding Ratio Conjecture can be written 
as 

(1-18) ^ Y = ^Y + 0(e"'^'’(-^)))), 

/6.F f&T 

for some <5 > 0. The exponent —5c{f) in the error term here indicates a power saving over the size 
of the main term. There are specihc instances in which more is known about this size of this error 
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term for averages of ratios of L-functions in various families - these are described in Section 15.21 
where there is a full discussion of this error term. 

In [23 Conrey and Snaith demonstrate how to extract information about the n-point correlation 
function of zeros from the n-over-n ratio of the Riemann zeta function. They use Cauchy’s theorem 
to express a sum over the heights of the zeros of the zeta function, 7 ^, as 


/(7l,---,7n) 

0<7i,...,7n<T' 


(1.19) 


C', 


/ • • • / ■ • • > TT ^( 1/2 + Zj) dzi... dzn, 

{2m)^ Jc Jc C 


where C is a positively oriented contour which encloses a subinterval of the imaginary axis from zero 
to T. The product of logarithmic derivatives of the zeta function can be obtained from a Ratios 
Conjecture by differentiating with respect to the a parameters in the numerator in an expression 
like (jl.lljl . This connection between ratios of L-functions and zero statistics is what we draw on 
in the current paper. We will not go further into the details of m as all steps of the method are 
covered explicitly in the present paper for matrices or L-functions with orthogonal or symplectic 
symmetry. 

We can see this brings us in a full circle back to the statistics of zeros of zeta and L-functions, 
as information about correlations of zeros are encoded in these Ratios Conjectures. We will see 
how this information is extracted in this present paper. First, in Section [2] we will use ratios of 
characteristic polynomials to derive a new form of the n-level density for the orthogonal group 
SO{2N). Section [3] follows, illustrating the same calculation for unitary symplectic matrices in 
USp{2N). Then the Ratios Conjectures of Conrey, Farmer and Zirnbauer are used in Sections [5] 
and [6] to conjecture the n-level densities for zeros of specific families of L-functions. In Section [7] we 
demonstrate why this new form of the random matrix theory n-level density simplifies immediately 
when the support of the Fourier transform of the test function (the equivalent of / in (jl.l9jl ) is 
restricted. Section [ 8 ] gives examples of 1- and 2-level densities explicitly. 


2. Eigenvalue Statistics of Orthogonal Matrices 

Let X he a 2N x 2N matrix with real entries X = (xjk)- We define the transpose matrix, 
X* = (xkj) and call X orthogonal if XX* = I. We let SO {2N) denote the group of all such 
orthogonal matrices with DetX = 1. All the eigenvalues of a matrix X G SO{2N) have absolute 
value 1 and can be written as , •'' > with 0 < 0 i, • • • , 0 ^ < tt. 

The result we are proving in this section is: 
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Theorem 2.1. 


( 2 . 1 ) 


The n-level density for eigenvalues of matrices from SO{2N) can be written as 
N 

E ,0,JdX 


L 


ji¥=jk'^i,k 


1 


(liri) 


E 


( 2 «) 


\M\ 


KULUM={1,--- ,n} 


7r\ n 


J*{—izK U izi) 


X f{zi,--- ,Zn)dzi ■■■dZn, 

where J*{—izK U izi) will he defined in \2.2(j\) and the sum in the right hand side involving K, L 
and M is over disjoint subsets of {1,... ,n}. 


2.1. Integral Theorem. In this section we will express the n-level density of eigenvalues of matri¬ 
ces from SO{2N) as contour integrals on the complex plane. We will see in the subsequent sections 
that moving the contours onto the imaginary axis (or, by a change of variables, the real axis) to 
arrive at Theorem 12.11 is non-trivial and will take some work. 


Theorem 2.2 (Integral theorem for SO{2N) matrices). Let C_ denote the path from —5 — ni up 
to —6 + TTi and let (7+ denote the path from 6 — 7ri up to d + iri. Let f be a 2TT-periodic, holomorphic 
function of n variables such that 

f (%i) • • • = f r • • , ±Gjn) 


( 2 . 2 ) 

Then 


(2.3) 




N 


/ (^ii) • • • ) %„) dXso{2N) 


(27ri) 


KULUM={1,--- ,n} 


j(jK JqLUM 


J {ZK u -zl) f {izi, • • • , iZn) dzi--- dz„ 


where zk = {zk '■ h E K} and —zl = {—zi ; I E L}, J^k Jqlum means we are integrating all the 
variables in zk along the C+ path and all others up the C- path and 


(2.4) J{A)= [ l[(-e-)^{e-ndXsoi2Ny 

JS0{2N) Ax 

Here Ax{e°‘) = det(/ — e^X*) is the characteristic polynomial of X; K,L,M are finite sets of 
distinct integers; A is a finite set of complex numbers; and d-^ 5 o( 2 X) indicates integration with 
respect to the Haar measure. 

Note that the sum over the eigenangles in (j2.3j) is not restricted to a sum over distinct indices. 
These extra terms, which do not figure in Theorem EU will cancel out the residues caused by 
moving the contours onto the imaginary axis in Section 12.41 
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X is an orthogonal matrix, so it has 2N eigenvalues , • ■ ■, such that 0 < 0i,..., < tt. 

The proof follows by considering g(z) = Ax(e^}, the characteristic polynomial of X. Using Cauchy’s 
theorem we can show that 

N t 

(2.5) ^ ,0jn) = -XAT / —{zi,---Xn)f{—,---—)dzi---dzn 

■ \zm) Jc-n g i % 

Jli’” )Jn — 

where C is a positively oriented rectangular contour around a interval of the imaginary axis of 
length 27r e.g. with corners ±z 7 r ± 5. By the periodicity of /, the horizontal integrals cancel and 
we use the functional equation and a change of variables from z ^ —z to get the result. 


Proof of Theorem\2fM The characteristic polynomial of X, 


(2.6) 

giz) = Ax (e^) 


N 

(2.7) 

= (1 - e"e*^^)(l - 


i=i 


has seros at a, = + 2nmi, m€Z,l<j<N. Then by Cauchy's Theorem (f) has poles 

at Zj = ^iOj + 2TTmi with residue 


( 2 . 8 ) 


9 (u) 

9 (zj) 



So let C be the positively-orientated rectangle with vertices ±(5 ± iri, where 5 is a small positive 
number. We can express the sum 


(2.9) 


( 2 . 10 ) 


N 

2” E /w..'" ."tJ 

jl,--- Jn = l 


N 

~ /(fl^il)’’’ Xndjn) 

in-jn=iee{i,-i}" 



The 2” term in front of the sum comes from counting all the possible combinations of +0* and —Oj. 
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We average this over X G SO{2N), and apply a change of variable from Zj —Zj. Let dG 
denote / {izi, ■ ■ ■ , izn) dzi, ■ ■ ■ dzn 


( 2 . 11 ) 


N 

2"/ E fiOnr--,0jJdX 

~ To W / ({^1 ) ■ ■ ■ ) ^n}) dG 

(Zm) J(jn 


where 


( 2 . 12 ) 


j(A)= f n 

JsO{2N) Ax 


By the periodicity of the function /, the horizontal segments of the contour cancel. Only the 
vertical paths need to be considered. 


(2.13) 





,Zn})dG 




n 

J {{zi,--- ,Zn})dG. 


So the expression is a sum of 2” terms, where each term is a n-fold integral each variable being 
integrated over C_ or C+. Another way to write this is 


(2.14) 


(2.15) 


/C’ 


J {{zi,--- ,Zn})dG 

_ ri 


dG 


KVJL={1,--- ,n} 


[ E (- 1 )'"' n / 


jeK' 


c+ Ax ^ ^ 


ic. 




where the sum over iL U L is a sum over disjoint sets. 
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For each variable Zj on C-, we replace —e (e “) with —2A^ + (e“) (obtained by 

differentiating the functional equation Kx{e^) = Kx{e~^)). 


(2.16) 


(2.17) 


(2.18) 




J ({zi, • • • ,Zn])dG 


c+ ^ 


7so(2N) 

={l,-,n} 

xnX 


,A' 


(e-^0 

c+ Ax ^ ^ 


f y: (- 1 )'""^' n [ 

JsO(2N) xuLUM 

^ (2A)l^l / [ J{{zkU-zl}) 


dGdX 


dG 


KULUM 


where the sums over A U F U M are sums over disjoint sets. 

By the definition of J{A) in (12.121) this gives us the required result. 


□ 


We can then use the Ratios Theorem to find alternative ways of expressing J (A). 


2.2. Ratios Theorem. Here we write down the average of ratios of characteristic polynomials of 
the even orthogonal group. The following form is rewritten into set notation from the result of 
Conrey, Forrester and Snaith [22] - they credit the original result to a preprint by Conrey, Farmer 
and Zirnbauer [20j . 

Take finite sets A and B of complex numbers. Let A > |R| and 5R(/3) > OV/3 G B and consider 


(2.19) 

( 2 . 20 ) 


R{A]B) 


I naeA^x(e 

JsO{2N) ll/3gsAx(e P) 

V / ^(^- U {A\D),D- U iA\D))Z{B,B)Y{B) 

V U {A\D),ByY{A\D)Y{D-) 







ORTHOGONAL AND SYMPLECTIC n-LEVEL DENSITIES 


13 


where D = {—a : a E D}, A\D = {a £ A, a ^ D} and 

(2.21) z{x) = 

1 — e ^ 

(2.22) y(^) = n z(2a), 

a&A 

(2.23) Z{A,B)=Wz{a + ^) 

osA 

P&B 


We can express J{A) as defined in (12.121) in terms of R{A] B) 


(2.24) 


aSA 


B=A 


So by differentiating the Ratios Theorem we can obtain a theorem about averages of logarithmic 
derivatives. 


Theorem 2.3. Let A he a finite set of complex numbers where 3f?(a) > 0 for a £ A and |A| < N, 
then J{A) = J*{A) where 


(2.25) 


(2.26) 


J{A) = / n i-e-^^ie-^dX 

SO{2N) 

^ -2N y; 5 

J*{A)=^e (-l)I^L 

DCA ' 

R 


I Z{D,D)Z{D-,D-)Y{D) 

Y{D-)Z^D-,D)^ 


llHniWr 


A\D=WiU-UWr r=l 
\Wr\<2 


where the sum over Wi is a sum over all distinct set partitions of A\D. 


(2.27) 


HoiW) = 


Ez (a—(5) 2 : (a+(5) 

z{a—S) z{afi-S) 


SeD 


{a + a) 


- ^(2a) W = {a} C A\D 

W = {a, a} C A\D 

W = 0 
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and 

=1 

1 — e ^ 

(2.28) 

(2.29) 

YiA) = n z(2a), 

aeA 

(2.30) 

Z{A,B)= llz{a + P). 

aeA 

h&B 


The dagger adds a restriction that a factor z{x) is omitted if its argument is zero. 

The main steps in the proof are pulling the differentiation inside the sum over subsets D C A, 
and then separating the differentiations in cases by whether a ^ D ov a ^ D. The differentiation 
by a G D are relatively straightforward, but to do the remaining differentiations we need to use 
logarithmic differentiation, as in the unitary case [26]. We do not include the details here as it is a 
simpler case of the proof of Theorem 15.21 presented later. 


2.3. Residue Lemma. In this section, we will locate the poles of J*{A) and calculate the residue 
of these poles. 

It is clear the only possible poles of J*{A) are when a = —j5 for some a, fi & A, or when a = 0. 
We need to know what the residues are at these poles. 

Theorem 2.4 (Residue Theorem for SO{2N) matrices). Let A be a finite set of complex numbers, 
let a*,/3* € A and K = and let J*{A) as defined in 112 . 26 \) . Then there is a simple pole 

when a* = —/3* and the poles cancel at a* =0. 

(2.31) Res^^ = J* (^A' U {,5*}) + J* (^A' U {-/S*}) + 2NJ* 

(2.32) Res (J*(kl)) = 0. 

a *=0 


We will consider the cases of (I2.31h and (I2.32h separately. First we will prove that the pole at 
a* = —j3* is simple and show that (|2.3ip holds. 

Since proofs of this structure are going to come up again, we define 

(2.33) X = A\{a*,^*} 

(2.34) D =Dr\X 

(2.35) {A\D)' = {A\D) n A 
where A\D = {a G A, a ^ D}. 

Definition 2.5. We say that the meromorphic functions Q and H of several variables have property 
Pf if the following four conditions hold for a continuous single variable function /. 

PI: If oT,j3* G A\D, then Q{D) is independent of a* and /3* and 

(2.36) H{D,W) = !. ^ ifW-{a,/3} 

\ 0(1) otherwise 
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(2.37) 


P2: li a* ^ D and j3* E A\D, then Q{D) is regular when a* = —/3* and 

^+ 0 ( 1 ) iiw = m 

0 (1) otherwise 


H{D,W) = 


(2.38) 


H{D,W) = \ 


P3: If a* E A\D and /3* E O, then Q{D) is regular when a* = —/3* and 

^ + 0(1) AW = {a*} 
otherwise 

P4: If a* E O and /3* E D, then Q{D) = + 0{l))Qi{D) where 

Qi{D) = QiD'){l - {a*+n{fin + H{D',{a*})\o^,=.0.+H{D',{f3*})) 

(2.39) +0(|a*+/3*p)) 
and 

H{D, W) = H{D', W) - (a* + /3*) W U {a*})a*=-/3* + 77(0', W U {/?*})) 

( 2 . 40 ) +0{\a* + P*\^). 

We will choose our f{x) dependant on what family of matrices we are considering. 

Lemma 2.6. If Q{D) and H{D,W) have property Pf and 

( 2 . 41 ) r{A) = ^ PDiA\D) 


DCA 


(2.42) 




DCA A\D=U, Wr 

IWrl<2 


(where Pd is defined by comparison with then 

(2.43) Res JJ*(7l)) = + J* U {/3*}) + J* U . 

Proof. If a*,j3* E A\D then 

(2.44) 


Pd{A\D) = Q{D) E n H{D, Wr 

r 

( 


A\D={_\^Wr r 
lWrl<2 


(2.45) 


(2.46) 


1 


{a* + 

1 

{a* + /3*fi 


Q{D) 




{A\D) =U, Wr 
\Wr\<2 


\ 


/ 


+ 0 ( 1 ) 


Pd{{A\D)) + 0{1), 


by PI. So Res (Pd(^\O))=0. 
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If a* G D and fi* G A\D, then 


(2.47) 

Pn{A\D) = Q{D) Y1 

X{H{D,Wr) 


A\D=U, 1 

/Ur r 



\Wr\<2 



(2.48) 

~ a* 

E 

J]i7(Zl,ITr)+0(l) 


(A\D)'=\J^Wr 

r 


a* P /3* ^ 

\Wr\<2 


(2.49) 

-P}iiA\D)' 

) + 0(1), 


by P2. Thus 


(2.50) ResJP d{A\D)) = P^,^^_^,^{{A\D)'). 

a.*^p* L j 


If a* G A\D and /3* G D, then, similarly to the previous case, 


(2.51) 


(2.52) 


(A\D)'=[J^Wr ^ 
\Wr\<2 


by P3. This allows us to conclude that Res {P£){A\D)) = P£){{A\D)'). 
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If a* G -D and (3* G D, then 


Pd{A\D) = Q{D) E n H{D, W, 

A\1 

r 

Q{D 


(2.53) 


(2.54) 


A\D=[J^Wr r 
\Wr\<2 


{a* + / 3*)2 




A\D=\J^Wr r 
\Wr\<2 


+ 


QjP') 

a* + ( 3 * 


E , Wr) fin+HiD',{n)+HiD',{-n) 


A\D=\J^Wr r 

\Wr\<2 




H{D', Wr U {/?*}) + II(D', Wr U {-/?*}) 


H{D',Wr 


;PD'i^\P) 


+ 0 ( 1 ). 

1 

[a* + /3*)2 

+ {finPD iA\D) + P^, iiA\D) U in) + Pd' HAD) U {-/?*})) 

+ 0 ( 1 ), 


by P4. Therefore 


(2.55) Res (Pd(A\0)) = finPn iA\D) + Pj,, ((^\0) U {/?*}) + Pd' HAD) U {-/?*}) • 

a*^p* 


We obtain the result of this Lemma by combining the results for the four cases. 


Res (J*(Gl)) 


(2.56) 


^ finPDiA\D) + P^, iiA\D) U in) + Pd' HAD) U {-/?*}) 

DCA 

a*,/3*eD 


+ E Piyui-f.) {a\d)') 

DCA 

a*eD,l3*eA\D 

+ Pj, ((Gl\0)') 

DCA 

l3*eD,a*eA\D 


(2.57) 

(2.58) 


= fin Y, PDiA\D)+ Y paad) 

DCA' DCA'u{/3*} 

+ ^ Pd{A\D) 

dca'u{-i3*} 

= fif3*)J* + J* [a' U {/?*}) + J* U {-/?*}) . 
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Note that the term in (12.541) features Pj^i{A\D) which is equal to Pd \^{A\D)j when 

a*, j3* G A\D, so the term cancels with the term from the first case, confirming that the 

pole at a* —>■ /3* is simple. 

□ 


It is straightforward to show that if f{l3*) = 2N, 


(2.59) 

and 


Q{D) = Y. 


-2N Y ^ 


Sen 




DCA 


IZ{D,D)Z{D-,D-)Y{D) 

Y{D-)Z^D-,Df 


(2.60) H{D, W) = Hd{W) = < 




(a + a) 


W = {a, a} C A\D 

hh = 0 


then PI to P4 hold, and so we have proved (|2.31l) . 

Now we will prove (12.321) . Again we define a general property. 

Definition 2.7. We say that the meromorphic functions Q and H of several variables have property 
R if the following two conditions hold. 


(2.61) 


Rl: If a* G A\D, then Q{D) is independent of a* and 

^ + 0(1) ifW = {a*} 
otherwise 


H(D,W) = { 


R2: If a* G D then 


(2.62) 


(2.63) 


and 


H{D,A\D)\o,,=^ = H{D ,A\D) 


Q{D) = ^Q{D^) + 0{1). 


Lemma 2.8. If Q{D) and H{D,W) satisfy these properties and 
(2.64) r{A) = Y Pd{A\D) 


(2.65) 


DCA 

DCA A\D=U^ Wr r 

IWrl<2 


(where Pd is defined by the above equation), then 


( 2 . 66 ) 


Res(J*(A)) = 0. 

a *^0 
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Proof. This is proved in a similar manner to Lemma 12.61 We have that if a* E A\D then by R1 


(2.67) 


R^s^(Pn(A\D)) = ^Pn (( 

A\D)') . 

If a* E D then by R2 



(2.68) 


Res (Pn(A\D)) = 

CK*—>-1) Z 

{A\D). 

Combining the two cases gives 

US 


(2.69) 

R«(J-(A))= ^ ipD((A\0)') 

- E Ipd'A\D) 



DCA 

DCA 



a*eA\D 

a*CD 

(2.70) 


= E ((■^\»)') - 

E \pd-a\o) 



d<za' 

DCA' 

(2.71) 


= 0. 



With the definitions (j2.59p and (j2.6n|l it is straightforward to show that properties R1 and R2 
hold, thus proving (I2.32|) . 


2.4. n-level Density of Orthogonal Matrices. So far we have expressed the n-level density of 
eigenangles of orthogonal matrices in terms of contour integrals. Recalling that C_ denotes the 
path from —6 — ni up to —6 + ni and C+ the path from 5 — vri up to <5 + vri, we take a / be a 
27r-periodic, holomorphic function of n variables such that 


(2.72) 

Then we have shown 


/ , • • • ,ejj = f , • • • , ±ejj 


L 


N 


(2.73) 


/ (^ii) • • • ) dXso{ 2 N) 


(27rz) 


A'ULUM={1,--- ,n} 


/ / J {zk'J-zl) f {izi,-■ ■ ,iZn)dzi-■ ■ dZn 

Jc¥ 


where zk = {zk : A: E iL} and —zl = {—zi : / E L} and /qlum means we are integrating all the 
variables in zx along the C_|_ path and all others down the C- path. The sum over K U L U M is 
over disjoint unions and J(A} defined as in (j2.4|] . 

We can then deduce the following. 
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Lemma 2.9. Let n < N 


(2.74) 


N 

2”/ E fiOnr-- ,e,JdX 

^ (2iV)l^l 


{2TTiy 


KULUM={1,--- ,n} 


/ / J* {ZK U -Zl) f (izi, • • • , iZn) dzi--- dZr, 

JcS 


Proof. We know that J{A) = J*{A) when |A| < N and 5R(a) > 0 Va G A by Theorem 12.31 This 
condition is clearly met by {zk U —zl). □ 

The next step is to move these contour integrals onto the imaginary axis. In order to do this, 
we first need some new notation. 

For given n, 0 < R < n, let 


(2.75) 



N 

E 

ji,;-- Jn=l 


For fixed sets K, L, M such that K Li LU M = {1, • • • ,n} let l m integral in Lemma 

12.91 with — i? of the integrals shifted onto the imaginary axis. All the integrals on the imaginary 
axis are principal value integrals. 


(2.76) 


Tn,R 

^f,K,L,M 




IC 


{LUM)n{l,-- 


■,R} 


r{zKU-ZL) 


X f{izi, • • • , izn)dzi ■ ■ ■ dzRdzR+i • • • dzn- 


We can express Lemma 12.91 in the new notation. 


(2.77) 


(27ri)’"2’ 


Is 


SO{2N) 




/(»,„■■■ ,e,„)dx 


A'ULUM={1,--- ,n} 


We will be using Sokhotski-Plemelj Theorem |84] [75], which states that we can shift a contour of 
integration onto line passing through a pole of the integrand. We interpret the resulting integral as 
a principal value integral and gain half the residue of the pole. For further details see any standard 
text on complex analysis, or a section on functions of a complex variable in a text such as [3]. 

We will now prove 
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Theorem 2.10. With the notation defined above, 0 < R < n 


(2.78) 


(27ri)-2" [ ,e,jdx 

JS0{2N) 

= E 

A'ULUM={1,--- ,n} 




Proof. We will prove this by induction, noting that we have already proved the case when R = n 
for all values of n. 

Firstly we will prove the base case when n = 1. We have already proven it is true for n = 1 and 
7? = 1, so we just need to show Equation (|2.78p holds for n = 1 and R = 0 


(2.79) 

(2.80) 
(2.81) 


(27ri)2 [ f {6,,)dX 

JS0{2N) 

= {2m)2 f 

JS0(2N) 


ISO{2N) 


K.L.M 


KVJLUM={1} 

[ j*{e)de- 

Jc+ 


'C- 


J*{-9)d9- 


2Nd9 


ic- 


By Theorem 12.41 the only pole of J*{9) is at 0 and has residue 0. 


(2.82) 

(2.83) 


/ r{9) + J*{-9) + 2Nd9 

J—in 

A'ULUM={1} 


We now move on to the inductive step. Assume Equation ()2.78l) is true for n = p — 1 and 
0 < R < p — 1 and consider the case when n = p, 0 < R < p. 

We will proceed by induction on R. We have already proved that Equation ()2.78p holds if R = p 
so take that as the base case. Assume that Equation (|2.78p holds if n = p, and R> S so that 


(2.84) 


{2Tri)P2P [ 

Js 


f{9,„--- ,9,^)dX 

SO{2N) 

Y (2iV)l^l/P;J, 

A'ULUM={1,--- ,p} 

Y (21V)l^l 

A'ULUM={1,--- ,p} 




(2.85) 


/ 27 r rin r p 

-i^'" J-in ■«> ^ 

X f{izi,--- ,izp)dzi,--- ,dzs,dzs+i---dzp. 
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Consider moving the zs contour on the right hand side of this equation, from just off the imag¬ 
inary axis onto the imaginary axis. It picks up residues from the variables whose contours have 
already been moved to the imaginary axis. Take t > s = S, then there are residues when zs = —zt 
and S,t G K oic S,t G L and also residues when zs = Zt t > S and S G L,t G K oic S G K,t G L. 

Consider fixed K,L,M. If S G K, then the residue of J*{zk U —ZL)f{izi, ■ ■ ■ ,izp) at zs = 
is 


( 2 . 86 ) 


ITT 


2NJ* ( 


Zr^' U -Z, 


■) + r( 


^K' ^L'u{t} 


+ J ( ^K'u{t} ^L' 


X ,izs-i,izt,izs+i,‘ 


where K = K {A — {S,t}),L = L (1 {A — {S', t}). The residue is multiplied by in rather than 
27ri because the zs contour is moving precisely onto the imaginary axis, so it only gives half the 
residue of a contour going completely around the pole by the Sokhotski-Plemelj Theorem |84] [75]. 
As / is symmetric in all its variables, the residue is the same when t G L and t G K as the residue 
is symmetric in t. 

The other residues are when S G L. Then zs appears in J*{zk U —ZL)f{izi, ■ ■ ■ ,izp) with an 
additional minus sign and is on the C- contour and hence being integrated around the pole in a 
clockwise direction. These two minus signs cancel leaving the same residue as in Equation (|2.86l) . 

Returning to our calculation in Equation (I2.85p . we note S is fixed and we split the sum into 
three different cases: S G K,S G L and S G M. We then consider the residues that result when 
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t > S and t G M as t varies from S + 1 to p. 


(2.87) 


,p} 


M 


KULUM={ 1 ,--- ,p} 

+ E E (2Ar)l^l 

t=S+l KULUM={ 1 ,- ,p} 


' —277 J C 


,s-l} y^(LUM)n{l,--- ,5-1} 


iTT {2NJ*{zj^i U —Zj^i) 


+ J [Zj^t U ^ ^L' 

X /(izi,--- ,izs-i,izt,izs+ir ■ ■ ,izp)dzi, ■ ■ ■ ,dzs-i,dzs+i ■ ■ ■ dzp 

+ E E (2iv)i^i 


t=S+l AULUAL={1,--- ,p} 
SeL,t^M 


' —in J C_ 


j^Kn{l, - .S-l} 


Kn{i, - .S-i} /^(LuM)n{i,--- ,S-i} 


—in (2NJ*{zp^i U —Zp^i ) 


+J {zj^i U z:J^l^J^^^) + J U Zj^i^ 

X f{izi,--- ,izs-i,izt,izs+i, - ■ ■ ,izp)dzi,--- ,dzs-i,dzs+i ■ ■ ■ dzp 


This expression can then be simplified to 

tP,S-1 
^f,K,L,M 




( 2 . 88 ) 


AULUM={1,--- ,p} 

+ 47rf E E (2Ar)l^l 

t=s+i k'ul'um={i,--- ,p}-{s,d 


• —ITT J C 


K ,5-1} l^{L uM)n 


[ / (2NJ*{Zr^lU—Zrl) 

J„(L UM)n{l, -- , 3 - 1 } '' ' ^ ^ ^ 


+ J {Zj^i U ■*" (%'u{d ^L' 

X f{izi,--- , zs-i, Zt, zs+i, ■ ■ ■ 'iZp)dzi, ■ ■ ■ ,dzs-i,dzs+i ■ ■ ■ dzp 


where the factor of four comes from the fact that there are 4 different combinations of S', t ^ M 
(i.e. S G L and t G K, S G L and t G L etc). All of the unions over sets are disjoint unions. 

Notice that f{zi, • • • , zs-i,zt, zg+i-, ■ ■ ■ ,zt,---Zp) has a repeated variable. It is one of the terms 
we want to remove from ED so that we have a sum over distinct indices. With this in mind we 
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relabel the variables zi, ■ ■ ■ , zs-i, zs+i, ■■■ Zp as zi, ■■ ■ , Zp-i. We define a function 

(2.89) gt,s(zi, ■ ■ ■ ,Zp-i) = f(zi, ■ ■ ■ ,zs-i,zt-i,zs, ■ ■ ■ ,zt-i,- ■ ■ ,Zp-i) 

(2.90) = f(zi, ■ ■ ■ ,zs-i,zt,zs+i, ■■■ ,zt,--- Zp). 

It is obvious that gt^s is the function / with zs = zt, for some t > S. 

We also note that for functions h of sets K, L and M 

h{K U {m},L, M) + h{K, L U {m}, M) + h{K, L,MU {m}) 

, , A'ULUAf={l,... 

(2.91) 

h{JC,C,M). 

KXJCyjM={l,-" ,m} 

This allows us to rewrite Equation ()2.88l) 

KULUM={1,— ,p} 

KULUM={1,--- ,p} 

t=S+l KULUM={1,--- ,p-l} 

By our inductive hypothesis on n, Equation (|2.85l) holds for n = p — 1 and R = S — 1. 

A'ULUM={l,...,p} 

= Y (2iV)l^l/^ 


(2.92) 


(2.93) 


tP,S-1 


KULUM={1,- ,p} 
P 

+ Am 

t=s+i 


i Y ( 2 ^*)^ ^ f Y^ ^ 9t,s{0hr-- ,0jp-i)dXso{2N)- 

^_C,1 JSO(2N) 


It is clear from the definition of in Equation (|2.75l) that 

V— •\TI,R ^ — •\TI.R —1 —^77.—1,/2—1 

(2.94) ^ ,9^)=Y ,en)+ Y 9t,R{du--- , 0 ^). 

t=R+l 

Considering now the left hand side of Equation (|2.85|) , we see that 
{27rir2P [ Y"’’' >''' , dj,)dXsoi2N) 

JS0{2N) 

{2mr2P [ /(0i, • • • , ep)dXso(2N) 

JS0(2N) 

(27ri)^2^ / 'Y Y^ 9t,s{9i,-■ ■ i9p-i)dXso{2N)- 

4 so f 2 AT) ,777 1 ^ 


(2.95) 


+ 


ISO{2N) 
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Equating Equations (I2.93p and (j2.95p . we see that 

( 2 . 96 ) A'uLu«=(..-,ri 

= (2nir2- f ,e,)dXso,m)- 

JS0{2N) 

Thus we have completed the inductive step in R by showing that Equation (I2.78P holds for 
n = p, R = S — 1 a it holds for n = p,R = S. As we have already shown Equation (|2.78p holds for 
n = p, R = p, we have now shown that it holds for all R when n = p. 

This completes the inductive step in n, as we have shown that Equation (I2.78h holds for all 
0 < R < n when n = p if it holds for all 0 < ii < n and n = p — 1. We already proved that it holds 
for all 0 < i? < 1 when n = 1, so by induction Equation (I2.78p is true for all n and 0 < R < n. □ 


Using Theorem 12.101 we can prove our main theorem on n-level density for orthogonal random 
matrices 


Theorem 2.11 (n-level Density of Orthogonal Matrices). 

/ ,9jjdx 

JS0{2N) 

1 


(2.97) 


{Iniy 


(2iV)l^l 

X ,Zn)dzi ■ ■ ■ dZn, 


J*{-iZK U iZL) 


where J*{—izK U izi) is defined in 112.26\) and the sum notation is defined at ^2.75\ l and indieates 
that the sum is over distinet indiees. 


Proof. We have proven that 


(2.98) 


(2.99) 


(27ri)^2" 


JsO{ 2 N) 


fio. 






A'ULUM={1,--- ,n} 


Y (2A^)l^l 

A'ULUM={1,--- ,n} 



n 

J*{zK'J-ZL)f{iZl,--- ,iZn) 


X 
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where the integral here is a principal value integral. Changing variable to 9i = izi etc 


( 2 . 100 ) 


( 2 . 101 ) 



J*{-izK U izL)f{zi, ■ ■ ■ 



Y, {2N)\^\r{-izKyJizL)f{zu--- 

KULUM={1,--- ,n} 


) Zn) 


) Zn)- 


Now 


( 2 . 102 ) 


^ (2iV)l^lj*( -izx'JizL) 

KULUM={1,— ,n} 

= Y {2N)\^\j*{iZKl^-iZL) 

KULUM={1,— ,n} 


and / T^jn) — f (=t%i, • • • , ±6*j„) so each integral from —tt to vr is double the integral from 

0 to TT. We change into this form to allow easier comparison with other correlation functions in 
other contexts. 


(2.103) 



) (^jn)dXso(2N) 


= 2 


n 



Y {2N)\^\r{-iZKUiZL) 

KULUM={1,— ,n} 


X f{zi,--- , Zn)dzi ■ ■ ■ dZn- 


All that is now required is to show that there are no poles on the path of integration. As / is 
holomorphic, we just need to check that 

(2.104) Y {2N)\^\r{-izK^izL) 

KULUM={1,--- ,n} 

has no poles at zi = Z 2 (this argument applies equally well to any other pair of z's by symmetry). 
We do not need to check for a pole at zi = —Z 2 as z* > 0, for 1 < i < n on the path of integration. 
A given J*{—izK U Izl) has a simple pole at zi = 2:2 if 1 £ A', 2 G L or 1 G L, 2 G AT. So 

Res I y {2N)\^\ J*{-izK U izi) I 

21 = 2:2 \ I 

\KULUM={l,-,n} ) 

(2-105) = y (2fV)l^l Rjs 

A:u£UA^={ 3,-,n} ^ 

+ d* (-^^a:u{^ 2} U izcuizr})'] 
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as ReSx=yf{x,y) = -Resx=yfi-x,-y). Therefore if (Y1kulum={i,-^ izl)^ 
has a singular set, it has complex dimension less than n — 1. However, this implies the singular set 
is trivial (see Corollary 7.3.2 in [62]). 

□ 


3. Eigenvalue Statistics of Symplectic Matrices 


We also want to look at families of L-functions with symplectic symmetry, so we will repeat the 
calculations we have done in Section [2] with USp{2N) . 

Let Y = (yjk) be a 2N x2N matrix and define the transpose matrix = (ykj)- T is a symplectic 
matrix if YZY^ = Z where 

^=(-1 ‘ o ) 

and /at is the identity matrix of size N. The group of symplectic matrices USp{2N) is a subgroup 
of U{2N) the group of unitary matrices, just like SO{2N), . All the eigenvalues of a matrix 
X ^ USp{2N) have absolute value 1 and can be written as , •'' > with 0 < 0i, • • • , < tt. 

The result we are proving in this section is: 

Theorem 3.1. The n-level density for eigenvalues of matrices from USp{2N) can be written as 


(3.2) 


L 


N 


7r\ n 


1 


(27ri)’' 


(2iV)l^l f / 1 '^USp{2N)i izK'ZizL) 
,rMA.f- \Jo / 


KULUM= 


X f{zi,--- ,Zn)dzi ■■■dz„ 


where defined at \3.12\) and the sum in the right hand side involving K, L 

and M is over disjoint subsets of {1,... ,n}. 


3.1. Integral Theorem. In this section we will proceed exactly analogously to the SO{2N) case 
and express the n-level density of eigenvalues of matrices from USp{2N) as contour integrals on 
the complex plane. We will perform exactly the same manipulations to move these contours onto 
the imaginary axis in order to arrive at Theorem 13.11 

Theorem 3.2 (Integral Theorem for Symplectic Matrices). Let C- denote the path from —b — rxi up 
to —b + ni and let C+ denote the path from b — iri up to b + iri. Let f be a 2'K-periodic, holomorphic 
function of n variables such that 


(3.3) 


/ = f (±6»j^,--- ,±6jJ. 
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Then 


L 


N 


(3.4) 


/) / (^Ji) • • • , ) dXusp(^2N) 

^ ^ (2iV)l^l 


{2'Kiy 


K\JL\JM={1,— ,n} 


/ / JuSp(2N) i^K U -Zl) / (^^1, • • • , 

JcS 


where zk = {-Sfc ^ k G iC} and —zl = {—^ G L}, J^lum means we are integrating all the 

variables in zk along the (7+ path and all others up the C- path and 


(3-5) JuSP(2N){^) — [ n '^)dXusp(2N)- 

JuSp( 2 N) Ax 

Here K, L, M are finite sets of integers and A is a finite set of complex numbers and dXusp(^ 2 N) 
indicates integration with respect to the Haar measure. 

The proof is identical to the proof for Theorem 12.21 


3.2. Ratios Theorem. Here we rewrite the theorem of Conrey, Forrester and Snaith [22] in set 
notation. Take finite sets A and B, where 7/ > |R| and consider 


(3.6) 


(3.7) 


where D 

(3.8) 

(3.9) 

(3.10) 


7?C/Sp(2Ar) (A; H) — / 

Ju 


{—a : a G D} and 


nagAAx(e °‘) 

USp(2N) 0/3e_B Axie-d) 

E 


dX 


-2iVE.,z.-5 Iy{A\D)Y{D) 


DCIA 


Y{B) 


IZ{D- U {A\D),D- U {A\D))Z{B,B) 
Z{D- U {A\D),Bfi 


z(x) = - - 

^ ^ 1 - e-^ 

Y{A) = n ^(2a) 

aeA 

Z{A,B)= llz{a + fi). 

aeA 

0eB 


Note that the main difference is the factor 
orthogonal case. 




y(B) 


Y{A\D)Y{D) 


in the 
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As in the orthogonal case, it is clear that Jusp( 2 N) = IlaeA ^-^t/S'p( 27 V) So the next step 

is differentiating Rjjgp(^ 2 N){^] B). 


Theorem 3.3. Let A be a finite set of complex numbers where 3ft(a) > 0 for a G A and |A| < N, 
then Jusp( 2 N){A) = Jhsp{ 2 N)i^) 


(3.11) 


(3.12) 


JuSp{2N){^) — j W_ '^)dXuSp{2N) 

TJ CteA 


USp(2N) 


J, 


USp{2N) o)= E 
DCA 


-2N E <5 
e (-1)' 

R 


IZ{D, D)Z{D-,D-)Y{D-) 
Y{D)Z^{D-,Dy 




A\D=WiU-UWRr=l 

\Wr\<2 


where the sum over A\D is a sum over all set partitions and 




/ z'{a-5) /(a+5)^ /,'2ai 

z(a+5) j z^) 

W = {a} C A\D 

(3.13) 

Hd{W) = < 

— [yj + ®) 

W = {a, a} C A\D 



1 

k. 

W = $ 

and 




(3.14) 


= 1 

1 — e ^ 


(3.15) 


y(A) = H z(2a) 




aeA 


(3.16) 


Z(A,B}= \[z{a + l3). 




oGA 

/3gB 



The dagger adds a restriction that a factor z{x) is omitted if its argument is zero. 
The proof is as in the orthogonal case. 


3.3. Residue Theorem. We need to know what the residues are at the poles of 

is clear that the only possible poles are when a* = —fi for some other /3 G A, or when a* = 0. 
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Theorem 3.4 (Residue Theorem for Symplectic Matrices). Let A be a finite set of complex num¬ 
bers, let a* ,j3* ^ A and X = A — {a*, j3*} and let Jijsp( 2 N)^X defined in Theorem \3.A . Then 

{jhsp[2N){A)) = JuSpi2NM' U m) + JuSp^2nM' U {-(3*}) 

+ ) 

5=0 {:^USp{2N){X) = 0 - 

Note that this residue formula is identical to that for J*{A) in the orthogonal case. 

Proof. The proof is very similar to that of Theorem 12.41 It follows by showing that the properties 
PI to P4 and R1 and R2 hold and then applying Lemma 12.61 and Lemma 12.81 □ 


(3.17) 

(3.18) 


3.4. n-level Density of Symplectic Matrices. Using Theorem 13.31 we can deduce 
Lemma 3.5. Let n < N. Then 


(3.19) 


^ fiOj,,--- ,6,JdX 


(27ri) 


A'ULUM={1,... ,n} 


f f 

Jc^ Jc 


J, 


qLUM 


USp{2N) 


(ZK u -zl) f (izi, • • • , iZn) dzi--- dZn 


Proof. We know that JuSp{2N){X — '^uSp{ 2 N)^X when |4.| < N and 5R(a) > 0 Va € ^ by Theorem 
ESI This condition is clearly met by {zk U —zl). □ 


Clearly as JijSp{ 2 N) same poles as J*{A) and an identical recursive formula for the 

residue, we can immediately deduce the following theorem. 


Theorem 3.6. Let 0 < R < n. Then 


(3.20) 


where 


(27ri)’^2" 


IUSp{2N) 


E ti.R 

m.. 




KULUM={1,--- ,n} 



n 

E 

in---jn=i 


(3.21) 
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and 


(3.22) 


I 


n,R 


f,K,L,M 


-r L 


,-R} 

+ 




(LUM)n{l,--- ,R} 


JuSp{2N)i^K^ Zl) 


X f{izi, • • • , iZn)dzi ■ ■ ■ dzRdzR+i ■ ■ ■ dzn- 


for fixed sets K, L, M such that the disjoint K U LU M = {1, • • • , n}. Note: this is the integral in 
Theorem ro with N — R of the integrals shifted onto the imaginary axis and all the integrals on 
the imaginary axis are principal value integrals. 


The proof for this in the orthogonai case rehes soieiy on the formuia for the residue, so this resuit 
foiiows immediateiy from the proof for Theorem 12.101 

We can now state the n-ievei density for sympiectic random matrices. 


Theorem 3.7 (n-ievei Density of Sympiectic Random Matrices). 

2” [ ,ejjdx 

JuSp{2N) 

1 


(3.23) 


{2TTiy 


(2iV)"^' ( r) '^USp{2N)i~'^^K U izi) 

lA/T— '^"'0 / 


A'ULUM= 

{l-.n} 


X fizi,--- ,Zn)dzi■■■dZn 


where JijSp{ 2 N)^~^^K U izi) is defined at \3.1^) and the sum notation is defined at 1^2. 75[ ) and 
indicates that the sum is over distinct indices. 


Proof. We foiiow exactiy the proof for the equivaient theorem in the orthogonai case - Theorem l2.11[ 
Taking the resuit of Theorem 13.61 we use a change of variabie to deduce 


(3.24) 


Noting that 


(27ri)-2- [ ,9jJdX 

JuSv(2N) 


>USp{2N) 

rTT \ n 


([ ^ {‘^N)^^^Ju.Sp{ 2 N)i-'^ZK^iZL)f{zi,--- ,Zn). 




(3.25) 


^ (2iV)WI JuSp{2N)i U izi) 

KULUM={1,— ,n} 

^ {2N)\^\jl,Sp^,^)iiZKiJ-iZL) 

A'UI,UM={1,--- ,n} 


and / (0j^,-- - ,0j^) = / (±0j^,--- ,±0j„) gives the required resuit. Foiiowing preciseiy the same 
steps as in the orthogonai case, it can be shown that there are no poies on the path of integration. □ 
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4. L-FUNCTIONS 


4.1. L-functions: a definition. We begin by defining what we mean by L- function. We will use 
a definition which is distinct from the “Selberg class”, but is conjectured to be equivalent to it m- 


Definition 4.1. Let s G C. An L-fnnction is a Dirichlet series 


(4.1) 


= Z] 

n=l 


(In 


with Un = 0{rf) for every e > 0 (Ramanujan Hypothesis [23]) subject to three conditions: it has an 
Euler product, an analytic continuation and a functional equation. We will define these properties 
similarly to the definitions used by Farmer in [35| and Conrey et al. in m- 


(1) The Euler product is a product over primes, equal to the Dirichlet series defining the L- 
function. Formally, for $R(s) > 1 we have 

W 

( 4 . 2 ) 

p i=i 

where the product is over the primes p, and each |7p,j| equals 1 or 0. tc is called the degree 
of the L-function. 

(2) The functional equation is an equation relating the L-function in one half of the complex 
plane to the other [58|. There exists e, such that |e| = 1, and a function T l(s) of the form 

W 

(4.3) ri(s) = P(s)Q*nr(|+/«.■)> 

i=i 

where Q > 0,‘Si{pj) > 0 and P is a polynomial whose only zeroes in Re{s) > 0 coincide 
with the poles of L{s) such that 

(4.4) iL{s) = VL{s)L{s) 
is entire and 

(4.5) ^l(s) = e^L(l - s). 

(3) The analytic continuation property means that L(s) continues to a meromorphic function 
of finite order with at most finitely many poles. 


It is worth noting that the Riemann Zeta function is an L-function with functional equation, 
analytic continuation and Euler product. 


4.2. Elliptic Curve L-functions. An elliptic curve, E, over a field iL is a non-singular curve of 
genus one, defined over K, with a iL-rational point, O ^ E. E has a algebraic group structure with 
unit element O. For more on elliptic curves see [87], [M] [82] . 

E can be embedded in as a cubic curve given by a global minimum Weierstrass eauation|89j 

(4.6) E : + aixy + a^y = + a2X^ -|- a^x -|- 
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with coefficients oi, • • • ,03 E K. Let E be an elliptic curve over Q, then the Weierstrass equation 
can be re-written as 

(4.7) E : + Ax + B 

with A, B G Zi and discriminant A = —{A^ + 27B‘^). For prime p, we can look at the curve 
y 2 _ ^3 - 1 - _j_ mod p. If p f A, then this is an elliptic curve over Fp and we say E has good 

reduction mod p. 

(4.8) Ep ■. = x^ + Ax + B mod p 

is a curve over the finite field Fp. If p\A, it is a prime of bad reduction. There are two types: mul¬ 
tiplicative reduction if E/Fp has a node, and additive reduction if E/Fp has a cusp. Multiplicative 
reduction is said to be split if the slopes of the tangent lines at the node are rational, and non-split 
otherwise. 


We define the conductor, M, of the elliptic curve as 


(4.9) 



K 

II 

where 


'0 

if E has good reduction at p, 

(4.10) 

fp = " 

1 

if E has multiplicative reduction at p. 


2 

if E has additive reduction ai p, p ^ 2, 3, 



^2 + 6p 

if E has additive reduction at p = 2 or 3. 


Here Sp depends on wild ramification in the action of the inertia group at p of Gal{^) on the Tate 
module Tp[E). [82| We will consider only prime M. 


Then for p of good reduction define the integer Up by 
(4.11) \E{Fp)\=p+l-ap, 


where |Fi(Fp)| is the number of points of Ep in a finite field of cardinality p, including the point at 
infinity. We define Op for primes of bad reduction by 


(4.12) 


dp — 


0 

< 1 
-1 


E has additive reduction at p, 

E has split multiplicative reduction at p, 

E has non-split multiplication reduction at p. 


We would like to use these Op to define an L-function in such a way that there is a L-function 
associated with each elliptic curve E. 


The Hasse-Weil L-function [52] of E/Q is defined as 

(4.13) u-B.*)=n(i-y) 'n(i 

p|A ^ ^ ^ pfA ^ 


A+pi-2^ 


-1 
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where A = + 27i?^). We define an for all n by expanding this into a Dirichlet series. 

OO 

(4.14) = 

n=l 

Hasse proved this was holomorphic for Re{s) > | by showing |ap| < 2y/p (Riemann Hypothesis 
for finite fields, proved by Hasse for elliptic curves in 1931 m)- The Hasse-Weil conjecture says 
that L{E, s) extends to a entire function. Eichler [33] and Shimura [81] proved that for elliptic 
curves E/Q that have a modular parametrization, L{E,s) can be extended to a entire function 
satisfying a functional equation 

(4.15) A{E, s) = u!e-^.{E, 2 — s) 
where wg = ±1 and 

(4.16) A(E,s)=(^^^ T{s)LiE,s). 

The Shimura-Taniyama-Weil conjecture says there is a cusp form whose L-function is equal to 
L{E,s), and implies the Hasse-Weil conjecture. Wiles [88] and Taylor [86] showed that this is true 
for all elliptic curves with square-free conductor. 


In 2001, Breuil, Conrad, Diamond and Taylor m then showed this was true for all elliptic curves 
over Q. 


Theorem 4.2 (Iwaniec and Kowalski [52]). Let E/Q be any elliptic curve of conductor M. There 
exists a primitive cusp form 

OO 

(4.17) f{z) = ^\{n)n^e{nz) G S 2 {To{N)) 

1 

such that A(n) = ^ giving 

(4.18) L{E,s + h = LE{s) = f2^. 

2 ^ n * 

From this it is clear that the Hasse-Weil L-function, L{E, s) meets all our conditions for an 
L-function in Definition 14.11 except the normalization of the symmetry line to 5R(s) = 


Le{s) = L{E,s + ^) gives a normalized L-function where the functional equation relates s to 
1 — s, giving 

(4.19) = Y: 




n(i 

p|A 


X{p) 


-1 


n 

pfA 


Kp) ^ J_ 

S 


-1 


(4.20) 




1 


F 
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with functional equation 

(4.21) Le(s) = 

where uje is +1 or —1 resulting in an odd or even functional equation for Le and 

(4.22) ^e(s) = r 0 - Le (s) = w(.E)^e( 1 - s) 

is entire. 


Clearly Le{s) meets all the conditions in Definition 14.11 for our definition of an L-function. We 
have demonstrated that each elliptic curve gives rise to a well-defined L-function. 

4.3. Dirichlet L-functions. We start by defining Dirichlet characters. 


Definition 4.3. [H] Let q € N. A Dirichlet character x to modulus g is a function y : Z — )• C 
such that 

(1) y (mn) = x {fn) X (^) for all m,n G Z] 

(2) X (q) has period q; 

(3) X (^) = 0 whenever (n, q) / 1; and 

(4) X(l) = l. 

A character is called primitive if it cannot be induced by any character of smaller modulus. For 
any Dirichlet character of modulus d, Xd (n), we can define a function 


(4.23) 

L (s, Xd) has a Euler product 

(4.24) 


Lis,Xd) = 


Xd{n) 


n=l 


ri'^ 


M»,x.)=n(i-^) 


-1 


and if Xd is primitive L (s, Xd) has an analytic continuation to the whole complex plane and has a 
function equation [33]. Setting 

s + a (Xd) 

I i I 

TT 


(4.25) 
where 

(4.26) 


^ (s,Xd) = i - 


L(s,Xd) 


a (Xd) = 


0 Xd(-l) = l 

xd(-i) = -i. 

Then the functional equation can be expressed as 

jp{Xd) 


(4.27) 


C(1 - s,Xd) = 


Ljt=iXd(t)e d 


i {s,Xd) ■ 
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Specifically note when Xd is real and d > 0, the functional equation is 
(4.28) L{l-s, Xd) = (^) ("’. 

Thus if we assume Xd is a real, primitive Dirichlet character L{s,Xd) meets all our requirements 
of an L-function in Definition 14.11 

4.4. Grand Riemann Hypothesis (GRH). We have shown that each of these L-functions asso¬ 
ciated to an elliptic curve or Dirichlet character has a functional equation with a line of symmetry, 
similar to the Riemann zeta function and its critical line. 

Conjecture 4.4 (Grand Riemann Hypothesis [52]). All non-trivial zeroes of L-functions L{s) in 
the critical strip 0 < 3ft(s) < 1 lie on the critical line 3ft(s) = 

In this work, we will assume the Grand Riemann Hypothesis holds and that L-functions arising 
from elliptic curves or Dirichlet characters have all of their nontrivial zeros on their lines of sym¬ 
metry i.e. the zeros have the form ^ iw with a; real.jS] 


4.5. Twisting L-functions into Families. 


4.5.1. Families of L-functions of Elliptic Curves. Starting with a primitive L-function generated 
by an elliptic curve, we can create a family by twisting with real quadratic Dirichlet characters. 
These characters can be described as follows. We first define the Legendre symbol for prime p, and 
all integers d 

-|-1 if p f d and = d mod p solvable, 

= < 0 p|d, 

— 1 if p f d and = d mod p not solvable. 

We then extend this to the Jacobi Symbol for positive odd integers m, 

k\ / r] 

■(si). 

where m = p^ ■ ■ ■ p^". 


(4.29) 


(4.30) 


d 
PJ L 


-) - {-) 

mjj Vpi/. 


The Kronecker symbol, is an extension of the Jacobi symbol for all integers d, n such that 

b\ fr 


(4.31) 


ab\ fa 


cd J 


K 


K \d/ K 


K \d/ K 


such that = (^) j for positive odd m, 


(4.32) 


d 


K 


1 

-1 


d > 0 
d < 0 
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and 

1 d odd, d = ±1 mod 8, 

(4.33) f-j = < 0 d even, 

^ —1 d odd, d = ±3 mod 8. 

The Kronecker symbol is a real quadratic Dirichlet character. From now on we will use Xd{n) = 

ii)K- 

We construct our family of L-functions by twisting with these Dirichlet characters. We start with 
a fixed primitive L-function Le{s) generated from an elliptic curves E as described in Theorem 14.21 

(4.34) - 


^Eis) = ^ ■ 




Using the quadratic Dirichlet characters we define a family parameterised by d. 

Kn)Xd{n) 


(4.35) 

(4.36) 

where 

(4.37) 


LE{s,xd) = y] ■ 


n* 


/ ^ _ Hp)Xd{p) '^Mip)Xdipf ^ ^ 


pa 


P- 


,2s 


^m{p) = 


1 if p t 

0 otherwise. 


Each L-function LE{s,Xd) is actually the L-function associated with another elliptic curve 
the twist of E by d [21] . Therefore it has a functional equation 


(4.38) 

(4.39) 


Le{s. X.) = xa(-M)^E (^) (1 -«) 

= Xd{-M)u;E'ilj{s, Xd)LE (1 - S, Xd) ■ 


and all the other attributes we expect of L-functions. Of particular interest is that the Grand Rie- 
mann Hypothesis in Conjecture 14.41 also applies to these twisted L-functions, so we are assuming 
that all the nontrivial zeros are on the Re{s) = i line. 


We will restrict our family of L-functions with the constraint Xdi—M)u}E = 1, where M is the 
conductor and cue the sign of the functional equation of the initial elliptic curve Le{s). This then 
gives us a family of L-functions with even functional equations: 


(4.40) E+ = {LE{s,Xd) : Xd (M) loe = +1; d a fundamental discriminant, d > 0}. 

It is expected that the distribution of the zeros of any one of these L-functions infinitely high on 
the critical line is like that of the eigenvalues of random unitary matrices. But rather than fixing 
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the L-function and looking at the distribution up the critical line, we can fix a height on the line 
and average through the family of E~^ with even functional equations. 


This way we can look at the distribution of the zeros on the critical line closest to the real axis 
such as the distribution of the ith zero 


(4.41) 








Kp 


G [a, b] 


l^xl 


where F is the family of L-functions with conductor Cf, and Kp is a, constant depending on the 
symmetry of the family. In the example above the conductor is the parameter d. Fx are the 
L-functions with Cf < X and t( is the ii\i zero of L-function / above the real axis. 


Katz and Sarnak conjectured that the zero statistics of twisted elliptic curve L-functions with 
even functional equation would behave like the eigenvalues of SO{2N) matrices and that zeros of 
families of L-functions with odd functional equation should behave like eigenvalues of SO{2N + 
1) |57]. For example, for the family F = F~^ we expect 

(4.42) lim F)[(i,b] = lim meas s A G SO{2N) : - — —^ G [a, 6] 

X^oo ’ ’ N—ioo [ TT 

where meas G SO{2N) : ^ 5]| is the distribution of the ith eigenvalue of a matrix A 

varying over SO{2N). 

4.5.2. Families of Dirichlet L-functions. Recall that we defined L{s,Xd) for some real, primitive 
Dirichlet character of modulus d, Xd (n) as 

( 4 . 43 ) = 

n=l 

We can then define a family of L-functions = {L(s, Xd) ■ d a fundamental discriminant, d > 0}, 
where d (the conductor of the L-function) is the ordering parameter of the family. 

Katz and Sarnak m conjecture that the zero statistics of D~^ behave like the eigenangles of 
USp{2N) matrices, citing evidence from Ozluk and Snyder [71] and Rubinstein |7S] to support this 
conjecture. 


5. Zero Statistics of Elliptic Curve L-functions 

In this section we will consider a family of quadratic twists (with even functional equation) of an 
elliptic curve L-function. Their zeros are conjectured to behave similarly to the eigenangle statistics 
of matrices from SO{2N) with Haar measure. We will apply the method of |27j to conjecture the 
n-level density of the zeros of these L-functions. We will see that the steps mirror the calculations 
for orthogonal matrices in Section [2j The main result to be derived in this section (stated formally 
in Theorem 15.101 and conditional on the Generalised Riemann Hypothesis for L-functions in this 
family and the Ratios Conjecture discussed later in this section) is the following form of the n-level 
density for zeros of elliptic curve L-functions in a family. For simplicity at this stage we denote the 










ORTHOGONAL AND SYMPLECTIC n-LEVEL DENSITIES 


39 


family by J- ^ a member of the family is labelled E and is the ith zero of that L-function above 
the real axis. We find the n-level density is: 


X] Z] fi'yjuer--,ljn,£) 

^ ' KULUM ^ 

,n} 

X f{zi, ■ ■ ■ ,Zn)dzi ■■■dZn + O , 

where J’^i—izx U izL,—izM) is defined in (jb.Q.'lji for the family of quadratic twists that we are 
interested in and the sum in the right hand side involving K, L and M is over disjoint subsets of 
{1,..., n}. We note how the structure of this n-level density mirrors that of the equivalent statistic 
for eigenvalues of matrices from SO{2N) given in Theorem 12.11 


5.1. Ratios of Elliptic Curve L-functions. Let A = {oi,-- - ,ax} and B = { 71 ,-•• ,7 q} be 
sets of complex numbers and e G {1, —1}^. In order to derive the n-level density of the zeros near 
the critical point s = we consider ratios of twisted L-functions with even functional equation 
averaged over the family E~^ defined at ()4.40p . All sums of the form 0 < d < X in this section 
include only fundamental discriminants. 


(5.2) 


-Re («!, • • • • • • )7q) 


= E 

0<d<X 

Xd(-M)a}E=+i 


Le + ai,Xd) ■■■ Le{\ + ax, Xd) 
Le {\+li,Xd) ■■■Le {\+lQ,Xd) 


We use the expression given for this quantity in [ 2 T] ; details of the derivation when K = Q = 1 can 
also be found in [48] . 

This gives us the following conjecture, from “Autocorrelations of ratios of L-functions” EU, which 
we have written in set notation as in the random matrix calculations. 

Conjecture 5.1. (Ratios of L-functions of Elliptic Curves) Let A, B he sets of complex numbers 
such that 


Im (a ), Im ( 7 ) <C 


Va G A 


(5.3) 


V 7 G R 
Ma e A,x ^ B 
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Then 


(5.4) 


Re{A,B) 


E 

0<d<X 

Xdi-M)u)E=+i 


OggA (2 Xd) 


= E E 

0<d<A DCA 
Xd(-M)uE=+l 


( VM\d\ 

27r 


~‘^^seD ^ 

n 


SeD 


r(i + 5) 


X Ye (A, B, D) Ae (A, B,D) + 0 


where the sum over d is over fundamental discriminants, C{s) is the classical Riemann zeta function 
and 


(5.5) 

(5.6) 

(5.7) 

(5.8) 

(5.9) 


Z^{A,B) 

YciA) 

U-{A) 

U+{A) 

V{A) 


C (1 + « + /3) 

aeA,l3eB 

nC(l + 2a) 

a&A 


n 

a£A 


n 

aSA 


A(P) 1 \ 

j^I/2+q; pl-\-2(Y J 



A(p) 

pl/2+Q: 


+ 


1 


P 


l+2a 



Hp) \ 

pl/2+a J 


For D C A we define D = {—5; S G D} and we can write 


(5.10) 

Ye{A,B,D) 

(5.11) 

Ae,i {A, B, D) 


Ae,2 {A, B, D) 

(5.12) 


(5.13) 

Ae,3 {A, B, D) 

(5.14) 

Ae {A,B,D) 


{{A\D) U D-, {A\D) U D-) {B, B) (B) 

Z^ {{A\D) U D-,Bf Y^ {A\D) Y^ {D-) 

1 


Ye{A,B,D) 

TT W U- (B) 


+ 


n 


U+{B) 

2U+ {{A\D)UD-) 

VjB) 

V{{A\D)UD-) 



Ae,i {A, B, D) Ae,2 {A, B, D) Ae,3 {A, B, D). 


5.2. Discussion of the Error Term. Little is known in general about the error term in the 
Ratios Conjectures, although we list below cases where some specific instances of the conjecture 
are better understood. It is clear that considerable error might arise in several steps of the ratios 
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“recipe” described in Section 11.21 which may be of a size similar to that of the main term. The 
conjecture is that these errors will cancel each other out and there are many situations where there 
is evidence for an error of the form O as stated in [2T] and reproduced in Conjecture 15.11 

Note that in (15.411 the average over the family (the sum over d) is not divided by the number of 
terms in the sum as was the case at (jl.181) . Also, for the family of quadratic twists of an elliptic 
curve L-function, the log conductor c(/) defined in Section [L2] is (asymptotically for large d) 21ogd 
for the L-function Le{^ + a, Xd)i so the general form of the error term given at (jl.181) would here 
correspond to for some <5 > 0. 

In several instances the error term in the Ratios Conjecture has been tested by assuming the 
Ratios Conjecture and deducing from it the 1-level density: a result like (15.1|) with n = 1. This is 
then compared with a calculation of the 1-level density using rigourous methods. These rigourous 
methods require a restriction on the support of the Fourier transform of the test function /. 

Miller rigourously calculated the 1-level density for the symplectic family of quadratic Dirichlet 
characters arising from even fundamental discriminants d < X. He found agreement with the 
1-level density computed using the Ratios Conjecture up to an error term of for test 

functions supported in (—1, l)|68j. 

Miller and Montague tested the accuracy of the Ratios Conjecture with an orthogonal family. 
They considered the 1-level density of families of cuspidal newforms of constant sign and confirmed 
the accuracy of the conjecture up to an error equivalent to when the test function has 

support in (—1,1). When the test function has support in (—2, 2) they can show a power savings 
error [70] . 

Goes, Jackson and Miller et ah showed similar results for the unitary family of all Dirichlet L- 
functions with prime conductor. They found agreement to the Ratios Conjecture predictions with 
a square-root error if the support of the Fourier transform of the test function is in (—1,1) and a 
power saving error for support up to (—2,2) [IT]. Similarly Fiorilli and Miller [36| calculated the 
1-level density for a family of Dirichlet L-functions and by explicitly computing lower order terms 
showed that the error term of type was the best possible. 

Huynh, Miller and Morrison tested the Ratios Conjecture prediction for the 1-level density of a 
family of quadratic twists of a fixed elliptic curve with prime conductor. They found agreement up 
to an error term of size A“ 2 ~ for test functions supported in (—<7, it)[19|. 

In the case of moments (ratios with no denominator), for the family of real, quadratic Dirichlet 
L-functions, a sequence of works [85] IH?)] IM] 190] indicate that in this case the third moment has 
a term of size 0{X‘^/^~^^ ) in place of the 0(A^/^+^) in (|5.4p . The size of this term was also 
extensively investigated numerically by Alderson and Rubinstein [2]. Very recently Florea m has 
found explicitly a term which in this notation would be equivalent to 0(A^/^+^) in the hrst moment 
of quadratic Dirichlet L-functions in the function field setting. 

Thus we have settings where there is support for the size of the error term being the square 
root of the size of the main term, as stated in the original Moment and Ratios Conjectures, such 
as (15.4p . and there are instances where the square root error term seems too optimistic. As we 
still lack a complete understanding of these terms, we will continue in this paper with the original 






42 


A.M. MASON AND N.C. SNAITH 


statement of the Ratios Conjectures with a O error term, but bearing in mind that there 

may be instances where this breaks down. 


5.3. Differentiating Ratios of Elliptic Cnrve L-functions. Using Conjecture 15.11 we want to 
find the logarithmic derivative 


Theorem 5.2. Assume Conjecture 15.11 Let A be a set of complex numbers such that 


(5.15) 


1 

logX 


<C 3? (a) < - 
4 


Im (a) <C ^ 


Mae A 
Mae A 


Then Je {A) 
(5.16) 


J%{A) + 0{X^I^+^), where 

Je (^) = n 

0<d<X aeA 
Xdi-M)u}E=+i 


J'E ih + 

Le (I + a,Xd) 


and 


(5.17) 


where 


j’Ea) = 


E E 

0<d<X DCA 

Xd{-M)ujE=+i 


i VM\d\ \ 
[ 2 - ) 


- E 25 


5GD 


n 

5eD 


rji-s) 

r(i + 5) 


^ Z^{D-,D-)Z^{D,D)Y^{D) 

y zl{D-,DfY^{D-) 

X (-1)1^1 Ae {D,D,D) 

R 

A\D=WiU-UWiir=l 


(5.18) 


(5.19) 


(5.20) 


Hd {Wr) 


Hd (Wr) 


AD,i (Wr) 


Hd (Wr) + (Wr) + Ad, 2 (Wr) + ^D,3 (Wr) 

YlseD (1 + a — <^) — y (1 + a + <5)^ — ^ (1 + 2a) 

< (1 + oi + 02 ) 

1 

0 

\ 


Wr = {a} 

Wr = { 01 , 02 } 
Wr = ^ 

|lUr| > 3 


n 

WSiWr 


A 

dw 


log AE,i (A,B,D) 


B=A 


and Ae-, Z(^,Yq etc defined as in Coniecture \5.1\ Ae(A,B,D) is an analytic function. The dagger 
adds a restriction that a factor (^(1 + x) is omitted if its argument is zero. The sum over d is over 
fundamental discriminants. 
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The proof of this theorem involves exactly the same arguments one would use to prove Theo¬ 
rem [231 although in Theorem 15.21 there is the added complication of the arithmetic terms. 

A necessary step in this proof is to evaluate Ae 

Lemma 5.3. Let A,B be sets and D, a subset of A. Then Ae {A,B,D)\^^j^ = Ae {D,D,D). 
This is an analytic function. 


Proof. Ae {A,B,D) is analytic as we can tell from it’s construction in Section [5.11 In particular, 
it must be analytic when A = B. However it is not immediately clear that Ae,i{A,B,D) is 
convergent when A=B due to the C (1) terms coming from Z {AfD U D~,B). In order to avoid 
problems of convergence, we will rewrite Ae,i {A, B, D) as a product over primes and consider the 
expansion of each prime separately. Recall 


(5.21) 
where 

(5.22) 

(5.23) 

(5.24) 


^E,i (^, B, D) = 


' Z{{A/D)VdD-,BfY{A/D)Y{D-) 

Z [{A/D) U D-, {A/D) U D-) Z {B, B,) Y {B) 


Z{A,B)= (^{1 + a + fj) 

y(A)= nc(l + 2a). 

a&A 


Now we will define 

(5.25) 

(5.26) Zp (A, B) = Zp(\ + a. + fi) 

a&A 

h&B 

(5.27) Tp (A) =l[zp{l + 2a) 

aeA 

noting that Zp{x) does not have a pole at 1 and that fdp Zp (A, R) = Z {A, B) and Y\_pYp{A) = 
Y (A). This allows us to rewrite Ae^i to get 


(5.28) 


^E,i,p (^5 B, D) = 


I Zpi{A/D)UD-,B)^YpiA/D)Yp{D-) 

Zp {{A/D) U D-, {A/D) U D-) Zp {B, B,) Yp {B) 


^E,i {A, R, H) = ]^ Ae,i,p (A, B, D). 


(5.29) 
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Then if we consider the value of Ae,i^p{A, B, D) when A = B 


(5.30) 

(5.31) 


Ae,i,p {A, B, — 


/ Zp{D-,DfYp{D-) 

V Zp{D-,D-)Zp{D,D)Yp{D) 
= ^E,i,p {D, D, D ). 


(5.32) 


Ae,i (-4, B, Ae,i,p {D, D, D). 


We will leave this as a product over primes and instead consider the second part of { A, B, D), 
Ae,2 


(5.33) 


(5.34) 

(5.35) 
where 

(5.36) 


(5.37) 


and finally 


(5.38) 


nrb 

p\M P 


1 ( U- (B) U+{B) 1 

1 + i \2U. {{A/D) U D~) 2U+ {{A/D) U D') p 


U+{B) 


1 ( U-jP) , U+jP) , 1 

l + l\2U_{D-) 2Up{D-) p 


-nrirUs 

p\M ^ P 

= Ae,2 {D, D, D) 


U-{A)=X{[l-^^+^ 


u+ (A) - n (^ + pTJ^ + ■ 


Ae,3{A,B,D)\p^^^= n 


V{B) 

V{{A/D)UD- 


- n 

p\M ^ ’ 

(5.40) = Ae,3{D,D,D) . 

where V {A) = O^eA (l “ • Putting these together, and remembering that Ae {A, B, D) is 

an analytic function, we get 


(5.41) 


Ae {A, B, D)\^^^ = Aep,p {D, D, D)Ae,2 {D, D, D) Ae,3 {D, D, D) 


( 5 . 42 ) 


= Ae {D,D,D). 
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With this lemma, we can then approach the proof of Theorem 15.21 The following steps are 
extremely similar to those necessary for the random matrix proof of Theorem 12.31 


Proof. Using Conjecture 15.11 we want to find the logarithmic derivative 


(5.43) 


(5.44) 


Je{A) 


E n 


0<d<X aeA 
Xd{-M)u}E=+l 


L'e (I + «.Xd) 

Le (^ + ot, Xd) 


d 




aSA 


B=A 


First we notice that if we substitute A in for B in Ye {A, B, D), we get 


(5.45) 


Ye{A,A,D) = 


IZ {D-, D-) Z {D, D)Y {D) 
Z{D-,DfY{D-) 


2 

The Z {D~,D) term in the denominator contains f (1), which is a pole. So Ye {A, A, D) is zero un¬ 
less D is the empty set as in the random matrices case. This means when we differentiate Re {A, B) 
by each a & D, all other terms are zero except the one where Z {D~, D) term is differentiated by 
every a G D. Let I3a represent the element of B that will be replaced by a when B is substituted 
for A. Expanding f and f about 1, it’s clear to see that 

^ (c(l + /3a-«)) 

Using this, we can see that 



(5.47) 


n 

a&D 


A 

da 


1 


Z{D-,Bd) 


Bd=D 


Zt {D-,D) 


where Bq represents the part of B that will be substituted for D, and Z^ means that we only 
include those terms of Z that are not equal to C (!)• All that remains now is to differentiate the 
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remainder by each ol G AjD. 


(5.48) 


where 

(5.49) 

(5.50) 

(5.51) 

(5.52) 

(5.53) 

(5.54) 


y. E 

0<d<A DCA \ 
Xd{-M)uJE = ^ 

X (-1)1^1 


( VM\d\\ 


2tt 


' E 25 

560 -|-rr(l-5) 

tiWTs) 


Z{D-,D-)Z{D,D)Y{D) 
Zt {D-,DfY{D-) 


n 

ya&A/D 


E 

da 


exp(// 


A,Bn 
D > 


\ 


Y (Bp) Z {Ba/d, Bpf Z {Ba/d, Ba/d) 
Z{D-,Ba/d)^ 


B=A 


+ 0 


H 


A,B 

D 


= log 


Z {A/D, D-f Z {A/D, A/D) 

\\ Z {A/D, Ba/dY Z {A/D, Bd)^ Y {A/D) 

X Ae {A,B,D)^ 

Y logC(l + a-7) + ^ E logC(l + a + /3) 

a&AID aeAJD 

X&D peAlD 

- Y logC(l + a + /3)-^ E logC(l + 2a) 


aeAlD 

/3eB 


a&AjD 


+ log {A, B,D)+ log Ae ,2 {A, B, D) 
+ log Ae ,3 {A,B,D) 


We now note that A H is a differentiable function of re G W, then 

/ \ 
d 


(5.55) 

where 

(5.56) 


n 


dw 


e^ = e^ 


Y H{Wi)... H{Wr) 

w=Wiu-uWr 




\w&A' 


The sum is over all set partitions of W into disjoint sets Wj. 
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We note that when we substitute D for Bu and A/D for then 

Z {A/D, D-f Z {A/D, A/D) Z {Ba/d, Bpf Z {Ba/d, Ba/d) Y (Bp) _ ^ 

Z {A/D, B^;nY Z {A/D, Bof , Z {D-,BA/DfY {A/D) 

and apply Lemma [531 to see that 


(5.57) 


(5.58) 

(5.59) 

(5.60) 


d 




yaeA/D 


\ 


Y {Bp) Z {Ba/d,Bd)^ Z {Ba/d, Ba/d) 
Z{D-,BA/Df 

R 


B=A 


Ae{D,D,D) E n HD{Wr), 

A/D=WiU-VJWRr=l 


where Hjy (VFr) is defined in Theorem 15.21 and plays the role of H{W) in (15.551) . This will allow us 
to simplify our solution after the differentiation. □ 


5.4. Residue Theorem for Elliptic Curve L-functions. Assuming the Riemann Hypothesis, 
the only possible poles of J%{A) are when a* = —/3* for some a*,j3* G A, or when a* = 0. We 
need to know what the residues are at these poles. 

Theorem 5.4. Let A be a finite set of complex numbers, let a*,13* G A and Al = A\{a*,/3*} and 
let J%{A) he as defined in Theorem \5.A Then 


(5.61) Res {J*e{A)) = J*e{A U {/3*}) + J*e{A U {-/?*}) - 

a*=—p* yj 


if f 1 


+ /3*) Xd ) Je{-^ ) 


where if {s,Xd) = Xd{-M)uiE{-^^"^ ' 

and 


comes from the functional equation for Le{s, Xd) 


(5.62) 


Res (J|;(A)) = 0. 

Q !*=0 


Proof. Let D C A, define Pd and Q as 

(5.63) J*e{A) = Pd{A\D) 

DCA 

R 

DCA A\D=W^\J-VJWRr=l 


(5.64) 
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where 


= E 


(5.65) 


0<d<X 

Xd{-M)a}E=+i 


f VM\d\ 

i 2-k 


- E 25 

5£D 


n 

<5eD 


r(i + j) 


IZ^{D-,D-)Z^{D,D)Y(;{D) 


V zl{D-,DfY^{D-) 

X (-1)1^1 Ae {D,D,D). 

He is defined in Theorem 15.21 and Ae, Z(^,Yq, etc, are defined as in Conjecture 15.11 

Using Q{D) in place of Q{D) and He{W) in place of H{D,W), conditions PI to P4 and R1 
and R2 (from Section 12.3p hold with 


(5.66) 


fin = - 


V'' A A o, /VM|dK r' 

^( 2 +^ -21og(^^)+- 


r' 


+ -(i-r) + ^(i+/3*). 


To show this we need to find expansions for Ae {Wr), He (Rr) and Ae (D) around a* = /3*. These 
lemmas will be given without proof, but complete proofs can be found in Amy Mason’s Thesis |65] . 

5.4.1. Taylor Expansions for Ae^x (Rr) j He (Rr) and Ae [D, D, D). Recall that 


(5.67) 


d 


AE,iiWr)= n -^^og Ae, i (A, B,D) 




for f = 1,2, 3. 


B=A 


Then we hnd 

Lemma 5.5. Let a*, (3* G H and D' = D/{a*, fi*}. For i = 1, 2,3 

d 

(5.68) 


da 


-AE,i (Wr 


a*=-l3* 


and 

(5.69) 


= i-AE,i (Wr + {/?*}) - Ae,x (Wr + {R}))L*=-/3* 


AD,xiWrn^_p,=A^,,iWr) 


A. B 

The second lemma we need for the Residue Theorem is to prepare for expanding 77^’ (Rr). 

Lemma 5.6. Let a*,(3* G D and D' = D/{a*, /3*}. For i = 1,2,3 

d 


(5.70) 

and 

(5.71) 


da 


-He (Rr 


a*=—j3* 


= i-HE{Wr + {n) 


- Rz5(Rr + K}))L*=-/3* 

HE{Wrn^_p,=H^, (Rr). 
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Most importantly, these lemmas show that when a*, /3* E D there is a Taylor expansion of 
(Wr) at a* = -/3* where (Wr) = Hd (Wr) + Ad,i (Wr) + Ad ,2 (Wr) + Ad ,3 (Wr) 


Lemma 5.7. Let Hd (Wr) etc. as defined in Equations \5.1^ fJ. 2(A Then 
ITd (Wr) = (Wr) 

(5.72) + (a* + n (-H^, (Wr + {fi*}) - (Wr + {a 

+ 0(|a*+/3*|2) 

Finally, we need the expansion of Ae (D,D,D) about a* = —fi* when a*, (3* El? 


a*=—^* 


Lemma 5.8. Let Ae {D, D) be as defined in Equations \5.11\ - \5.14\ and He {Wr) etc. as defined 
in Equations E73- \h.2(K 


III 


(5.73) 


Ae (D,D,D) = Ae(D ,D ,D 

i-(a*+n(^HD({n)-HD(m) 
+7{^({-/3*})-HD({-i 
+ 0(\a*+/3*\^) 


With these lemmas in place, we are ready to show that PI to P4 from Definition 12.51 hold for 
Hd, Q and / as defined above in Equation 15.661 

PI: Suppose a*,/3* E A/D, then Q(D), is independent of a* and j3*. 

As a* fi* the only pole in HD(Wr) comes when Wr = {a *so when 


(5.74) 


HD({a*,fi*})= (l + a*+n- 


Now (1 + x) 

— + O (1) and Ad (Wr) is analytic for all Wr- So 

(5.75) 


— {tt^^+ 0(1) W = {a*,fi*} 

Hd (Wr) = < ) 

I 0(1) otherwise. 


Hence the requirements of PI hold. 

P2: Now we consider the case when a* G D,fi* E A/D. Then Q (D) is independent of (5* and 
regular as a* — fi*. The only pole in Hd iWr) comes when Wr = {/3*}. Then 


Hd ({r}) = ^ f (1 + /3* - 5) - ^ (1 + r + <5) - y (2r). 

SeD 


c 


c' 




(5.76) 
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When J = a*, — ^(/3* + <5) is a simple pole with a residue of 1. (Wr) has no poles so 

(5.77) + 

lUtlj otherwise. 

This satisfies the conditions of P2. 

P3: Now we consider the case where a* E A/D, (3* E D. This is similar to the previous case, 
Q (D) is regular as a* = —/3* and the only pole in Hd/Wt) is when Wr = {a*}- 


(5.78) 


HoiWr) = 1 “*+/^* ^ 

10(1) otherwise. 


This shows that Q (D) and Hjy (Wr) satisfy P3. 

P4: We need to consider when a*,j3* E D. Recall that D' = D/{a*,j5*}. We already have 
expansions for Hd {Wr), Aj:)(Wr) and Ae{D,D,D) =: Ae{D) in Lemmas 15.5115.81 We want to 

find an expansion for Q {D) around a* = (3*. 


Q{D) := (-1) 


\ D \ 


(5.79) 


( VM\d\\ 

[ ) 


-2E, 


'5GD ^ 


n 

S£D 


r(i-<5) 

r(i + 5) 


I Z^{D-,D-)Z^iD,D)Y^iD) 

ZI{D-,D)^Y^{D-) 


= Q D 


/\ fVM\d\\ 


-2(a*+/3*) 


r(i-a*)r(i-/3*) 


(5.80) 


\ 2-71 I r(l + a*)r(l + /3*) 

c (1 + 2 a*) C (1 + 2 p*) C (1 + g* + /9*) C (1 - «* - 13*) 

C (1 + a* — f3*) C (1 — a* + (3*) 

C (1 - g* - <5) C (1 + W + ^) C (1 - /3* - <5) C (1 + /3* + (^) 
C (1 - g* + <5) C (1 + g* - 5) C (1 + /3* - <5) C (1 - /3* + '^) 


<5eD 


= Q D 


-1 


(g*+/3*)^ 12 


-- + 0(|g*+/3*|) 


X l-{a* + /3*) 


(5.81) 


Hj,, {{f3*}) + H^, ({-/3*}) 




+ 0(|g* +/3*p) 


Combining this with the results for He (Wr), AE{Wr) and Ae{D, D, D) = Ae{D). 
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(5.82) 


(5.83) 


(5.84) 


Q{D) = Q{D)Ae (D) 

= q(d 


^ + o {\a* + ^* 


{a* + /3 
l-(a*+r) 
il^' f 1 


12 


({r}) + ^n' ({-/ 3 *}) 


V' V2 


+ /3*, Xd 


+ 0{\a*+f3*\^) 


X Ae{D ) (^1 + (a* + n [-He' ({r}) + H^, ({/3*}) 

({-r}) + Hj,, ({-r})] + O {\a* + / 31 ' 
20^) 


-1 


(a* + p 
QiD' 


{a* + P* 

+ 0 ( 1 ) 




Combining this with Lemma 15.71 the conditions of P4 are met. Hence by Lemma 12.61 we have 
proved the first part of the Residue Theorem 15.41 

For the second part, we need to check that Q{D) and HE{Wr) meet the conditions R1 and R2 
from Definition 12.71 

Rl: Suppose a* E H/Zl,then Q{D) is independent of a*. The only pole comes from Wr = {a*}, 
where 


(5.85) 


^d(K)) = E7(i + 


d&D 


C C 

a* - d) - + a* + d) - + 2a*). 


The residue as x —0 of ^ (1 + x) is 1. So 


(5.86) 


HD{Wr) = 


Thus the condition Rl is met. 


^ + 0(1) 1F = K} 
0(1) otherwise 
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R2: Now we consider the case where a* E D. Then D), Ae (D) and Ae iWr) are regular 

at a* = 0 ,and He{A/D)\ a*=o = Hj^i{A/D). 


Q{D) = Q{D’) 


( 27r 


(5.87) 


\^\d\) 
C{l + 2a*)ll 


\- 2 «* r(l-a*) 


deD 


r (1 + a*) 

C{1 + a* + d)C,{l - a* - d) 
C(1 + a* — d)C(l + d — a*) 


The only pole comes from the C(1 + 2 q:*), which has a residue of This gives 
(5.88) q^ = ^5^) + 0(1) 

as required. Hence we can apply Lemma 12.81 and conclude the proof for the second half of the 
Residue Lemma [531 


□ 


5.5. n-level Density of Zeros of Families of L-functions of Elliptic Curves. In this section 
we will express the n-level density of zeros of families of L-functions of elliptic curves as contour 
integrals on the complex plane. We will define the n-level density fnnction, 5®(/), for the zeros 
for a family of twisted L-functions with even functional equation. 

(5-89) S^{f)= /(7ti,<i,--- ,7tn.d) 

0<d<X tl,--- ,tn>0 

Xd{-M)u}E=+^ 

where 7 *^^ is the height of the zero of LE{s,Xd) on the critical line, above the real axis. Note 
that the snm over the zeros in S^{f) is not restricted to a sum over distinct indices at this stage. 

Theorem 5.9. Let C- denote the path from —6 — ooi up to —<5-1- ooi and let (7+ denote the path 
from 6 — ooi up to 6 + ooi, where 6 is a small positive number. Let f be a holomorphic funetion of 
n variables sueh that 


(5.90) 

Then 


(5.91) 


/ r- - ,XjJ = f i±Xj ^, • • • ,±XjJ 


2 " (27ri)"5f(/) = 


/ / ( — 1 )'^'«7£; (zx U —Zl, Zm) 

X / (izi, • • • ,iZn) dzi--- dZn 


where zk = {zk '■ k E K}, —zl = {—zi : I E L}, J^k Jqlum means we are integrating all the vari¬ 
ables in Zk along the C+ path and all others along the C- path and 


(5.92) 


je(a,b)= y. 


0<d<X aeA 
Xdi-M)u}E=+i 


P&B 


ijj ^2 
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where ip (s, Xd) = Xd{-M)ujE ( 7 ^ 
are finite sets of eomplex numbers. 


\2s-l ■p(^-s) 

' j T{s+^) ■ ^ finite sets of integers and A, B 

The sum over d is over fundamental discriminants. 


This follows immediately from Cauchy’s Residue Theorem, in a similar manner to the proof of 
Theorem 12.21 

Using Theorem 15.21 and Theorem 15.41 we can replace our Je{A,B) with Jf.{A,B) + 
where we define 


rE{A,B)= Y. 

0<d<X 

Xd{-M)i^E=+^ 

TT r(l-^) I Z^{D-,D-)Z^{D,D)Y^{D) 
(5.93) Air(l + 5)y zl{D-,DfY^{D-) 

X (-1)1^1 Ae {D,D,D) 

R 

^ E U^DiWr) 

A\D= r=l 
WiU---UWji 


n jil+^^xd) E 

/3eB ^ DCA V 


/ VM\d\ \ 

V 2- i 


- 

5£D 


where He etc are as described in Theorem l5.2l Our previous definition of JeA) is the equivalent of 
JfiAj 0) in our new notation. All our theorems using the original definition have trivial extensions 
to encompass this new notation. 

Note, however, that this additional factor means that the residues found here may be more 

complicated than those in the random matrices case. In particular we need to consider whether 
/ 

the ^ + fi*,Xd) factors add additional poles, or affect the residues at the current ones when 

contemplating the multiple integrals in what follows. Recall that 

(5.94) I (1 + r. X.) = 2 log (^) 4 (1 - 4 (1 + 

Clearly the constant term is not a problem, so we only need to focus on the digamma terms. These 
have poles only at 0 , — 1 , —2 etc and are meromorphic elsewhere[T] so these are simple to deal with 
as long as |5R(/3*)| < We are already assuming this in Equation (I5.15P so it turns out the 
additional factor is harmless. 


We now have that 


(5.95) 


(2vrz)"2"5f (/) 

= E /k /l E Je{zk'J-zl,zm) 

KULUM= 0<d<X 

{1,•••,"•} Xd{-M)u}E=^ 

f {izi, • • • , iZn) dzi--- dZn + o 
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Note that we have pulled the error term outside of the integral as the size of the error is unaffected 
by integration with respect to Zi. 

Here we have assumed that / decays fast enough that the integrals are negligible outside the 
range given for the imaginary parts of the parameters in the Ratios Conjecture, Conjecture 15.11 
See the one-level density section of [26] for more details, but since it is not at all clear over what 
range the Ratios Conjectures hold, we do not pursue this further here. 

In line with the calculations with orthogonal matrices in Section [231 we will move all the contours 
one by one onto the imaginary axis. We recall from Section 15.41 that if a, /3 G A 


(5.96) 


Res J'^{A,B) 

a^—/3 


= JUA U {/3}, B) + JUA U {-/3}, R) - JUA , B U {/3}) 

where A' = ^\{q;,/3}. 

We define, for given n, 0 < R < n, 

rvD 

, , — ^n,R 

(5.97) 


E 


E 


jl,- Jn = l 


For fixed sets K, L, M such that the disjoint union K U L Li M = {1, • • • , n} let ^ be the 

integral in Equation (I5.95|) . excluding the error term, with — R of the integrals shifted onto the 
imaginary axis. All the integrals on the imaginary axis are principal value integrals. 


(5.98) 


jn,R 

^f,K,L,M 


I —ioo J C 


[ 

Kn{l, -,R} Jf- 


^{LUM)n{l,--- ,R} 


Je{zk U -zl,zm) 


X f{izi, • • • , iZn)dzi ■ ■ ■ dZRClZR+i ■ ■ ■ dZn 
We can express Equation (j5.95p in the new notation. 


(27ri)"2" 


E E”'”/( 7 : 




) ljn,d) 


(5.99) 


We will now prove 


0<d<X 

Xd{-M)u}E=+i 




KULUM={1,--- ,n} 


Theorem 5.10. Assume Conjecture 15.11 With the notation defined above, 0 < R < n, and d 
summing over fundamental discriminants, 


(27ri)"2" 




(5.100) 


0<d<X 

Xd{-M)uJE=+i 




A'ULUM={1,--- ,n} 


He 
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Proof. This proof follows identical steps to that of Theorem 12.101 so we will refrain from repeating 
all the explanation and instead just reference the relevant point in that proof of the SO{2N) case. 


The first step, in analogy to the steps leading up to (|2.83p . is to confirm the case n = 1, i? = 0: 


(5.101) 


(271*) 2 


0<d<X 

Xd{-M)LJE=+^ 



-f 

JE{z,9)f{-iz)dz + [ 

(5.102) 

Jc+ 

Jc- 

f 


— 

/ JE((I),z)f(-iz)dz 



Jc- 

(5.103) 

= E 



l—izjaz 




KULUM={1} 


We now consider the inductive step. Assume Equation (j5.100p is true for n = p—1 and 0 < R < 
p — 1 and consider the case when n = p, 0 < R < p. We will proceed by induction on R. We have 
already proved that Equation (I5.100p holds if i? = p so we take that as the base case. Assume that 
Equation (j5.100p holds if n = p, and R> S so that 

P,s 


{27ri)P2P Y, ,7ip.d) 


0<d<X 

Xd{-M)u}E=+i 


(5.104) 


AULUM={1,--- ,p} 




+ 0 ( 


AULUM={1,--- ,p} 
rioo rioo 


/ i<X) pioo r r 

-ioo " ' J-ioo - ^ 

X f{izi, ■ ■ ■ ,izp)dzi ■ ■ ■ dzsdzs+i ■ ■ ■ dzp + O . 


The residue structure is identical to the case in Theorem l2.101 Thus the only new information we 
need for this case is that for fixed sets K, L, M such that the disjoint union A'ULUM = {1, • • • ,p}, 
if S' G AT, then the residue of J%{zk U —zl, ZM)f{izi, • • • , izp) at zs = ±zt is 


(5.105) 


777 [-Je {^k' U -V’^MUp}) + Je 
+ J U —Zp^i , Z]Yj ^ 

xf{izi,--- ,izs-i,izt,izs+i, - ■ ■ ,izp) 

wherein = n (A — {S, t}), L = L D {A — {S,t}). 
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Thus, in analogy to equation (|2.88p . we find that the main term of (j5.104jl is equivalent to 


(5.106) 


A'ULUM={l,---,p} 

+ dvri ^ ^ (—1)1^1 

t=S+l K'yjL'vjM= 


r*2oo rioo 


'—200 J—iooJC, 


X n{i,- - ,5-1} 
+ 


X f{izi, • • • ,zs-i,zt,zs+i, • • • izp)dzi ■ ■ ■ dzs-idzs+i ■ ■ ■ dzp. 


All of the unions over sets are disjoint unions. 

Defining g as at (I2.90p . we can rewrite Equation (I5.in6p 




KULUM= 
{I,--- ,p} 


(5.107) 




M 


KULUM= 
{I.--- ,p} 


E E 


-1,5-1 
K,L,M' 


1=5+1 A'ULUM= 
{I,... ,p-l} 


By our inductive hypothesis on n, Equation (|5.104l) holds for n = p — 1 and R = S — 1. 




A'ULUM= 

{ 1 ,■■■,?} 


(5.108) 




5-1 

K,L,M 


KULUM= 


+ 47ri {2TTi)P ^2P ^ Y E^ \t,s{ln,dr-- ,ljp-ud)- 


1=5+1 


0<d<X 

Xtl( —= 1 
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The equivalent of (I2.95p results in the left hand side of Equation (15.104^ being written as 


{2mf 2P f i^h,d^ • • • > U.d) 


0<d<X 

Xd{-M)ojE=+l 


{2Tri)P2P ^ filjud, • • • , ljp,d) 


(5.109) 


0<d<X 

Xd{-M)ujE=l 


+ ( 2 «)» 2 » 


0<d<X t=S+l 

Xtj(—A^)t<2B=+l 


Equating Equations (I5.108p and (I5.109h . we see that 

AUI,UM={1,-" ,p} 


(5.110) 


{27ri)P2P ^ /(7ii.d, • • • , 7ip,d)- 


0<d<X 

Xd{-M)ujE=+i 


Thus we have completed the inductive step in R by showing that Equation (I5.100p holds for 
n = p, R = S — 1 if it holds for n = p, i? = S. As we have already shown Equation (|5.100l) holds 
for n = p, R = p, we have now shown that it holds for all R when n = p. 


This completes the inductive step in n, as we have shown that Equation (j5.100[) holds for all 
0 < R < n when n = p if it holds for all 0 < i? < n and n = p — 1. We already proved that it holds 
for all 0 < ii < 1 when n = 1, so by induction Equation (I5.100p is true for all n and 0 < R < n. □ 


We can state our hnal result for the n-level density of zeros of families of twisted elliptic curve 
L-functions. 

Theorem 5.11 (n-level Density of Zeros of Families of Twisted Elliptic Curve L-functions). As¬ 
sume Conjecture I5.il and the Riemann Hypothesis. Then 

^n,0 


E V-^ «i-, u 

2^ /(7il,d,--- ,7in,rf) 


0<d<A 

Xd{-M)u}E=+^ 


(5.111) 


(277) 


( roo\ n 

/) 

JO / 


KULUM 


J%{-iZK u iZL, -izm) 


X fizi,--- , Zn)dzi ■ ■ ■ dZn + O (^2+'"^ , 


where J'^i—izK U izL,—izM) is defined in L5. 9 A) and the sum notation is defined at ([g.751 ) and 
implies the sum is over zeros with distinct indices. The sum over d is over fundamental discrimi¬ 
nants. 
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Proof. This follows easily from the previous theorem. We have proven that 

^n ,0 


(27rz)”2" Y ’ 


(5.112) 


0 <d<A' 

Xdi-M)u)E=+l 


(5.113) 


E + o (x5+-) 

A'ULUM={l,---,n} 

E (-1)'"' 

A'ULUM={1,--- ,n} 


(^j ^ JE{zK'J-ZL,ZM)f{iZl,--- ,iZn)dzi---dZn + 0 1^X2+^'^ 


where the integral here is a principal value integral. 
(5.114) 


(2,rir2” E ^ 

0<d<X 

Xd{-M)u)E=+i 

Y (-1)'^' f / ) JE{-iZK'JiZL,-iZM)fizi,--- ,Zn)dzi---dzn + 0 

■iiniAd \J-oo/ ^ 


= I 


KULUM 

={l,-,n} 


(5.115) 


= ) Y i-'^)^^^JEi-'>'ZK'JiZL,-iZM)fizi,--- ,Zn)dzi---dZn + 0 (X 2 +A . 

\J — OO/ JV'I I r I 


KULUM 

={l,-,n} 


Now 


(5.116) 


Y {-lMU*E{-iZK^iZL,-iZM) 

KULUM={1,--- ,n} 


A:uLUAf={l,--- ,n} 

due to ^i^ + a,Xd) = ^i^-Oi,Xd)- In addition, f (9 ,9jU = f (±9j^,-" ,±%„)- Thus each 
integral from —oo to oo is double the integral from 0 to oo. 

(27)”2” ^ 

0<d<X 

Xd(-A4')t<2B=+l 
roo\ ^ 

Y (-lMUf;(-iZK UiZL,-iZM) 

KULUM={1,--- ,n} 


(5.117) 


= 2 " 


X f{zi, ■ ■ ■ ,Zn)dzi ■■■dZn + 0 . 
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All that is now required is to show that there are no poles on the path of integration. As / is 
holomorphic, we just need to check that 

(5.118) ^ (-1)1^1 UizL,-izM) 

A'ULUM={1,--- ,n} 

has no poles at zi = Z 2 (this argument applies equally well to any other pair of ^;’s by symmetry). 
We do not need to check for a pole at zi = —Z 2 as Zj > 0, for 0 < i < n on the path of integration. 
A given J%{—izK U iz^, —izm) has a pole at 21 = Z 2 if 1 £ AT, 2 G L or 1 G L, 2 G i^. So 


Res 

Zl=Z2 


(5.119) 


(5.120) 


(-1)1^1 U iZL, -izM) 

^A'ULUM={1,--- ,n} 

(-1)1^1 Res^ {JU-izKu{zi} U izLu{z2}^ -izm) 

KULUM={'i,--- ,n} 

+ >^e(“*^A'U{ 22} *^LU{zi}) —^Zm)) 


= 0 


as Resx=yf{x,y) = —Resx=yf{—x,—y). The same argument holds for all poles at Zi = ±Zj for all 
0 ^ ^ n. Therefore there are no poles on the path of integration, for the same reason as in 

Section 12.41 □ 


6. Zero Statistics of Quadratic Dirichlet L-functions 

Let Xd be a real primitive Dirichlet character. For any Dirichlet character of modulus d, Xd (n), 
we can define an L-function 

(«) = 

n=l 

L {s, Xd) has a Euler product 

( 6 - 2 ) = 

If Xd is primitive L{s,Xd) has an analytic continuation to the whole complex plane and has a 
functional equation [U]. 

When Xd is real and d > 0, the functional equation is 
(6.3) L{1-s,Xd) = {^) " (^Wd) • 

We can further define a family of L-functions D+ = {L(s, Xd) ■ d > 0 a fundamental discriminant}. 
All sums over 0 < d < A in this section include only fundamental discriminants. We expect that 
the n-level density of the zeros of the L-functions in D^, averaged over the family will correspond 
to that of the eigenangles of USp{2N). 
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Our plan is therefore to follow the recipe of Conrey, Farmer and Zirnbauer [2T] to express ratios 
of Dirichlet L-functions in appropriate terms and then mimic the calculations from Section [3] to 
determine the n-level density. 


6.1. Ratios of Quadratic Dirichlet L-functions. The result from “Autocorrelation of ratios of 
L-functions” EU, can be rewritten into set notation. 

Conjecture 6.1 (Ratios of Quadratic Dirichlet L-functions). Let A = {oi,--- ,a_fc} and B = 
{bi) • • • )7 q} such that 


(6.4) 


< 3f?(ai) < ^ 

logX < 4 

Ve > 0, Im (at), Im (jj) <C 


l<i<K 


l<j<Q 

I <i < K,1 < j < Q 


then 


(6.5) 


nfc=l (2 ^k-,Xd) 

0<d<X n^=l (2 + IqiXd) 

-E/6f/ r 

X Ydl (41, R, F) Adl (^, B,F) + 0 




0<d<X FcA 


where d is summed over fundamental discriminants and Ydl and Adl are defined for D C A as 
follows 


( 6 . 6 ) 


(6.7) 


I Z^ {{A\D) U D-, {A\D) U D-) Y^ {{A\D) U P-) {B, B) 

V Zd{A\D)UD-,BfYdB) 


Adl (41, B, D) = Ydl (4l, B, D)-^ J] ^ 

p p 


A V-{B) 

\2V- ((A\D) UD-) 


1 V+jB) 1 \ 

2V+{{A\D)UD-) pj 
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and 


(6.8) 

14(A)- n (l+ i+J 

aeA \ / 

(6.9) 

V-iA)=Yl{^- L) 

aeA \ / 

(6.10) 

D- = {-6:6 G D} 

(6.11) 

Z(^{A, B) = C (1 + ® + ^)) 

asA 

beB 

(6.12) 

y4A) = nC(l + 2a). 

asA 


6.2. Differentiating Ratios of Quadratic Dirichlet L-functions. Following the calculation in 
Section Ea we now wish to calculate 


(6.13) 

where 

(6.14) 


L (\ 


0<C.d^X 


jdl{a)= n X (9 ^ ^ n 7^^^^ 


d 


i^eA 


da 


B=A 


Rdl{A,B)= ^ 


YlkL{h + ^k,xd) 


0<d<X Hq -^ ( 2 + Xd) 

Theorem 6.2. Assume Conjecture \6.1\. Let A he a finite set of complex numbers where 

(a) < ^ yae A 

(6.15) log^ 4 

Ve > 0,lm{a) ,<€. X^~'^ ^a G A. 

Then we have Jdl{A) = Jf)j^{A) + O where 


(6.16) 


Jdl (A) = ^ 


0<d<X qSA 


L fl 

T 


+ Oi,Xd 


n^rfy 


0<d<X DcA 


SeD 


IZ(;{D-,D-)Z^iD,D)Y^{D-) 


(6.17) 


ZI{D-,D)Y^{D) 


X (-1)1^1 Adl {D,D,D) 
R 

A\D=WiU-UWr r=l 
\Wi\<2 


X 
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Hd (Wr) = Hd (Wr) + Ad {Wr 


(^( 1 + 7-/3)-7(1 + 7+ / 3 )) + 7(1 + 27) 1 +. = {7} 


/3eD 


C 


C 


( 7 ) (1 + 71+72) 


Wr = {71,72} 
Wr = ^ 


d 


(6.20) 71,, (1+,)= H —Adl{A,B,D) 


a&Wr 


B=A 


Adl{A, B, D) are defined in Conjecture \6.1\ and the sum over 0 < d < X includes only 
fundamental discriminants. 

Note that the conditions on A come from the conditions on Conjecture 16.11 

This calculation follows exactly the same lines as similar calculations and arguments used in the 
proof of Theorem 15.21 For full details, see the thesis of Amy Mason 


6.3. Residue Theorem for Quadratic Dirichlet L-functions. In this section, we will locate 
the poles of Jfij^{A) and calculate the residue of these poles. 

Theorem 6.3. Let A he a finite set of complex numbers, let G A and A' = A\{o, 6 } and let 

Jf)L{W) be as defined in Conjecture \6.A Then 


( 6 . 21 ) 


Res {JhLm = JhdA U (ri) + JhdA U {-fi*}) 

a*=—p* 




DL 


fi’DL V2 


1 


+ fi* Xd) Jdl{^ ) 


^ ^ r(—) 

where i^dl {s,Xd) = ( 7 ) ^ comes from the functional equation. 

^ \~T~' 


( 6 . 22 ) 


Res (J1,,.(A)) = 0. 

Q :*=0 


Proof. Let D C A, define Pd and Q as 

(6.23) Jh^A) = Pd{A\D) 


DCA 


(6.24) 


E ^ E n 


DCA 


A\D=\J^Wr r 
\Wr\<2 


(6.25) 
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where 


(6.26) 


< 3 ('D)= E E 

0 <d<X DcA 



E, 


Sen ‘ 


n 

Sen 


r(i-l) 

r(i + |) 


IZ^{D-,D-)Z^{D,D)Y^{D-) 


ZI{D-,D)Y^{D) 


X {D,D,D) 


and Hd is defined in Theorem 16.21 and Ae, Z^, Y^, etc, are defined as in Conjecture 16.11 

Replacing Q{D) with Q{D) and H{D, W) with Hd{W) in Section iTdl conditions PI to P4 and 
R1 and R2 hold with 

(6.27) = + 

This proof is similar to that in the elliptic curve L-function case and we will not reproduce them 
here. Full details can be found in Amy Mason’s thesis[65]. □ 


6.4. n-level Density of Zeros of Families of Quadratic Dirichlet L-functions. We will 
define the n-level density function, S^^{f) for the zeros for the family of L-functions D'*' = 
{L{s, Xd) : d > 0; d a fundamental discriminant}. 

(6.28) S^\f)= E .lu,d) 

0 <d<X tl,... ,tn>0 

where is the zero of L (s, Xd) on the critical line, above the real axis. Note that the sum 
over the zeros in S^^{f) is not restricted to a sum over distinct indices at this stage. 

Theorem 6.4. Let C- denote the path from —6 — ooi up to —6 + ooi and let C+ denote the path 
from 5 — ooi up to 6 + ooi, where 5 is a small positive number. Let f be a holomorphic function of 
n variables such that 

(6.29) /(d,i,--- ,d,J = /(±d,„... ,±%J. 


(6.30) 


! r^K I/^LUM 
KULUM='^^+ 


(_1)|A^| jDL u 


Then 

2-(2vri)"5f^(/) = 

E 

JLUA 

X / {izi ,... ,iZn) dzi... dZn 

where zk = {zk '■ h G K}, —zl = {—zi : / G Lj and J^k Jqlum means we are integrating all the 
variables in zk along the C+ path and all others along the C- path and define 


(6.31) 


jdl{a,b)= n ^( F +/ 3 wd ) 




/3gB 


IpDL 2 
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where il^DL{s,Xd) = (5) ^ with this definition being the obvious extension of Jdl{^) 


defined in Equation (|6.13l) . Here K,L,M are finite, disjoint sets of integers and A,B are finite 
sets of complex numbers. The sum over d is over fundamental discriminants. 


The proof of this trivially follows the same arguments as the proof of Theorem 15.91 Note that 
again we have an identity giving Jj:)l{A,B) = Jfi^{A,B) + O where A is a finite set of 

complex numbers where <C 3 ?(q;) < \ and Im{a) <C for a G A via Theorem 16.21 Again, 

all preceding results translate immediately to Jdl{A,B) from the original Jdl{A). Here 




-E, 


5gD ' 


0<d<X l3eB ^ ^ '' DcA 

r(i-l) jZ^{D-,D-)Z^{D,D)Y^{D-) 


(6.32) 


ZI{D-,D)Y^{D) 


^ n r (1 _g 5) 

X (-1)1^1 Adl {D,D,D) 

R 

X E 

A/D= r=l 
WiU---UWji 

with all terms as defined in Theorem 16.21 

The only possible poles are, as in previous cases, when a, /3 G A and a = 0 or a = —fi with 
residues 


(6.33) 


(6.34) 


B)) — 0 , 

Q ;*=0 


Res {JhiiAB)) = JhL U' U {r},R) 

a*=—p* \ / 

+ JhL (a' u {-/?*}, r) 
-JhL {A',Bu{fi*}) 


^ _ 1 P(’£') 

where fioL {s, Xd) = (f) ^ comes from the functional equation. Reusing our earlier notation 

defined in Equation (|5.97|) 

(6.35) 


E 


n.R 


E 


il,- ,jn = l 

jiJjXhk>R 


For fixed sets K, L, M such that KULUM = {1, • • • , n} let l m integral in Theorem l6.4l 

with — R of the integrals shifted onto the imaginary axis. All the integrals on the imaginary axis 
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are principal value integrals. 


(6.36) 


jU^R 


/ icxD rioo r 

-ioQ J —ioo J C 


,R} 

+ 




(LUM)n{l,--- ,R} 


Jdl{zk U -zl,zm) 


X f{izi, • • • , izn)dzi • • • dzRdzR+i ■ ■ ■ dzn- 
Then it can be proved, using the same ideas as the proof of Theorem 15.101 that 

Theorem 6.5. Assume Conjecture Id. 11 With the notation defined above, Q < R < n, and d 
summed over fundamental discriminants, 


(6.37) 


0<d<X 

KULUM={1,— ,n} 


This immediately gives us the result for the n-level density of the zeros of = {L{s, Xd) : d > 0} 
averaged over the family. We can state our final result for the n-level density of zeros of families 
Dirichlet L-functions with real quadratic characters. 

Theorem 6.6. (n-level Density of Zeros of Families of Quadratie Dirichlet L-functions) Assume 
Conjecture Id. 11 Then 

fhjudr-- ,7jn,d) 

0<d<X 

^79^ X] (-1)'^' f / ) J*DL{-iZKCliZL,-iZM) 

X f{zi, • • • ,Zn)dzi ■■■dZn + O 

with Jfi]^{A,B) as defined in \6.32\) and the sum notation defined at 116. 751) indicating a sum over 
zeros with distinet indiees. The sum in d is over fundamental diseriminants. 


Using the same argument as in Section 15.51 it is clear there are no poles on the paths on inte¬ 
gration. 


7 . n-LEVEL Densities with Restricted Support 

Currently it is extremely difficult to verify that number theoretic results on n-level densities 
agree with theorems from random matrix theory. These calculations always involve a test function 
(for example the function / in (15.111 ) which is sampled at positions of the zeros, and theorems can 
only be proven if the Fourier transform of this test function has its support restricted to a fixed 
interval around the origin. Rudnick and Sarnak [80] and Rubinstein |79] both use ad hoc methods 
to compare their calculations with support restricted to (—1,1), for L-functions with unitary and 
symplectic symmetry respectively. Gao m extended Rubinstein’s calculation to (—2,2), but was 
only able to verify his result for n < 3. Levinson and Miller |64j extended this to n < 7 and 
Entin et al.[34] verified it for all n, by-passing the complicated combinatorics that made previous 
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calculations so difficult and using a function field analogue. Still it remains clear that there is no 
consistent method of verifying these connections and that in every case the task is non-trivial. 

Conrey and Snaith |28j used their random matrix calculations from the earlier paper m to allow 
a more accessible n-level density comparison between random matrix theory and the Riemann zeros. 
They demonstrate that the way they write the random matrix theory n-point correlation formulae 
in m allows immediate identification of which terms remain when restrictions are placed on the 
support of the Fourier transform of the test function and so the simplified expression can be used 
to compare with a number theory calculation where the support is restricted. 

We give a sketch of the way in which the same ideas can be applied to the n-level density of 
eigenangles of orthogonal (see Section 17.Ih or symplectic matrices (in Section I7.2jl when similar 
restrictions are placed on the support of the test function. These sections give expressions for 
n-level densities of eigenvalues of orthogonal and symplectic matrices when the Fourier transform 
of the test function has any range of snpport, reducing the number of terms under consideration 
when matching with number theoretic results. These results should allow for easier and more 
straightforward verification of whether number theoretic results limit to random matrix expressions 
in the orthogonal and symplectic cases. Sections 17.11 and 17.21 describe the principle behind the 
calculations, while the details would depend on the particular family of L-functions and the set-up 
of the Fourier transform and any smoothing or weighting functions nsed; this is the subject of 
future work. 


7.1. Orthogonal. Recall from Section 12.41 that Lemma 12.91 tells us that for n < N and / a 
27r-periodic, holomorphic function of n variables such that 


(7.1) 

then 


(7.2) 


/ (Xi, • • • ,Xn) = f (±Xi, • • • , ±Xn) 


N 


2"/ E ,ejJdX 

^ (2iV)l^l 


(Itti) 


KULUM={1,--- ,n} 


/ / J* {ZK U -zl) f {izi, ■■■ , iZn) dzi--- dZn 

Jc^ 


with J*{A) etc as described in Section [231 


For consistency we rewrite this in terms of the notation used by Conrey and Snaith [28]. Firstly 
we redefine our idea of eigenangles so that we have an infinite set of points. For X a 2N x 2N 
orthogonal matrix with eigenvalues n = 1,... ,N, we will consider the statistics of the set of 

points 


(7.3) <■■■ <9o<0i<-- - <0 r< eR+^ < ■■■ 

where 


(7.4) 


dr+2kN — + 2TTk. 
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Take -F(xi, • • • , Xn) a holomorphic function of n variables such that 
(7.5) F (xi, • • • ,Xn) = F (±xi, • • • , ±Xn) . 

Let 5 > 0 and indicate integration from S — ioo to 5 + ioo and assume F{xi, ■ ■ ■ ,Xn) decays 
rapidly in each variable in horizontal strips. Then a simple extension of Lemma 12.91 gives 

CX> 

^ ,ejjdx 


L 


(7.6) 


SO(2N) 
1 


Jl,--- Jn = -0O 


(27ri) 


(2iV)l^l 


A:uLUM={ 1,... ,n} 


/ / / J* {ZK^ -Zl) F {izi,-■ ■ ,iZn)dzi-■ ■ dZn- 

7(5)1^! 7(-5)RI 


where 


(7.7) 


^ -2iV E d in, 
J*{A)=Y^ ''en 

DCA 

R 


'Z{D, D)Z{-D, -D)Y{D) 
y(-L»)Zt(F-,F)2 




A/D=VEiU-UlVflr=l 

\Wr\<2 

where HoiWr)-, Y, Z are defined in Equations 12.27112.301 

We apply this with F(xi,--- ,Xn) = rewrite / in terms of its Fourier 

transform, <l>j-, so that we can investigate the effect of limiting the support of 

' iN zi iNz„ 


f 


(7.8) 

where 

(7.9) 


27r 


’ 27r 


( 6 ) • • • ,^n)e 


2tt 


iN Zn^n 

27r 


d^i ■ ■ ■ d^ri 


e -- 


27r 


iN Zji^', 

271 


^ ' = ^NzAl + - + Nz„^r, 


Note that |3f?(2:i)| < 6 for all i and suppose that (^i, • • • , ^„) = 0 if |^i| + • • • + > 2q — e for 

some e > 0. That is, / has restricted support in the sense that arises in n-level density calculations 
in a number theoretical context. Therefore 


(7.10) 


^Nzi^l-\ - \-Nzni„ 


< e 


N5{2q-e) 


Now consider sets D where |I1| > q. Then, as |3ft(d)| = <5, Vd G D, the exponential term in J*{A) 
above is bounded by 


(7.11) 


3 2AJ]^g_pd 


< e 


-2N5q 
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Clearly then the product of these two factors in Equations 17.101 and 17.111 is also bounded 


(7.12) 


- \-NZn^n —2N Xlde D 


< e 


-NSt 


0 as 5 —>■ oo. 


It is significant that (I7.12p tends to zero as <5 goes to infinity because we know from Section 12.31 
that we can shift the (5) contours far to the right and the (—5) contours far to the left without 
encountering any poles of the integrand J*{zk U —zl) F{izi,... ,izn) of (17.61) . Thus if F decays 
fast enough as any z variable strays far from the real axis then (17.6h doesn’t change as the contours 
are shifted to infinity. It can be shown as in [28] that all factors in the integrand other than those 
considered in (I7.12h can be bounded and in this way we see that any term with \D\ > q in (17.7h 
contributes nothing to the integral. With this in mind we define 


(7.13) 


= E 


-2Af E <5 

e (-l)l'^l, 


DCA 

\D\<q 


' Z{D,D)Z{-D,-D)Y{D) 
Y{-D)Z^{D-,Df 


R 

^ E WHniWr). 

A/D=WiVJ-UWr r=l 
\Wr\<2 


When the support of / is limited to 2q, we can rewrite the result in Lemma 12.91 to give 

Theorem 7.1. Let 5 > 0 and indicate integration from 6 — ioo to 6 + zoo. Take F{xi, • • • , Xn) 
a holomorphic function of n variables such that 

(7.14) F (xi, • • • ,Xn) = F (±xi, • • • , ±Xn) . 


and assume F{xi, • • • , Xn) = / ' ■> decays rapidly in each variable in horizontal strips. 

With defined as in Equation \7.8[ let (^i, • • • ,^„) = 0 if |^i| + • • • + |^n| > 2q — e for some 
e > 0. Then 


(7.15) 


L 


SO{2N) 

1 


£ F{ej„--- ,ejjdx 


n,- ,Jn = -00 


(27rz) 


( 2 Ar)l^l 


A'ULUM={1,--- ,n} 


/ / / J* {ZK^ -zl)F {izi,--- ,izn)dzi---dz„ 

J(S)^ J{-5Y Jio)^ ^ 


where H£){Wr),Y, Z are defined in Equations \2.27\ - r2.d(J\ and Jq{A) defined in Equation \7.1‘J\ above. 


7.1.1. Special Case q = 1. We particularly want to consider the case g = 1 as the terms sim¬ 
plify dramatically, and should allow for for immediate identification with the terms survive the 
restrictions in the number theory case when the support is in (—2, 2). 
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Consider Theorem O in the case that <? = 1. Then in the sum within Jq{A) the only valid term 
is that which arises from D = 0 . 


(7.16) 


(7.17) 


where 


JKA) = Y, 


IZ{D,D)Zi-D,-D)Y{D) 


(- 1 )' 


DQA 

| D |<1 


Y{-D)Z\{D-,DY 


R 


Y XlHoiWr 


A/D=Wi\J-VJWiir=l 

\Wr\<2 

R 

Y 

A=WiU-UWii r=l 
\Wr\<2 


(7.18) 


770 (Wr) = 


( 2 a) Wr = {a} 

( 4 ) (« + / 3 ) Wr = {a,^}. 


7.2. Symplectic. We can calculate a similar result for the eigenangles of symplectic matrices. For 
X an 2N x 2N symplectic matrix, we extend its eigenvalues as before 

(7.19) ■■■e_R< 0^+1 < • • • < 00 < < • • • < < 0R+1 < • • • 

where 

(7.20) 0r.+2fcA = 0r + 27r7. 


Take F{xi, • • • , Xn) a holomorphic function of n variables such that 
(7.21) F (xi, • • • ,Xn) = F (±xi, • • • , ±Xn) . 


Let 5 > 0 and indicate integration from 8 — ioo to 5 + ioo and assume F{xi, ■ ■ ■ ,x„) decays 
rapidly in each variable in horizontal strips. Using Lemma 13.51 and the notation of Conrey and 
Snaith [28] we can see that 


(7.22) 


'USp{2N) 


Ji,- 


Y ,0jJdX 

,jn = -00 


1 

(2^ 


Y ( 2 iV)l^l 

A'ULUAL={1,--- ,n} 




duSp{2N) 


{zK'J-ZL)F{izi,--- 


,iZn)dzi ■■■dZn- 
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where 


(7.23) 


and 


r , ,(.11 - ^ j'ici / z(P.P)z(-p,-D)r(p) 

JvSp(2N)W- { 1) y y{-D)ZHD-,D)^ 


R 


Y. UHoiWr 


A/D=WiU-UWRr=l 

\Wr\<2 


(7.24) 


HniWr) = 


{ — {a — d) — — (a + d) 


deD 


(2a) 


j (« + /3) 


1 


Wr = {a} 

Wr = {a,/3} 

Wr = 0 


We apply this with F(xi, ■ ■ ■ , Xn) = / ^'^r) as in the orthogonal case we write / 

in terms of its Fourier transform. 

Suppose that (^i, • • • ,^„) = 0 if |^i| + • • • |^n| > 2g — e for some e > 0. Hence 

iNziCi 


(7.25) 


e -- 


\ ^ N5{2q-e) 


2tt 27r 

and considering sets D where |i7| > g in a similar manner to the orthogonal case above, we can see 


the exponential term in J^gp(^ 2 N) bounded by 


(7.26) 


3 


< e 


-2N5q 


Hence the product of these two factors in Equations 17.251 and 17.261 is bounded by 


(7.27) 

(7.28) 
Define 


(7.29) 


e -- 


27r 


iN Zn^r, 

2-k 


g— 2 A J2d£D ^ 


< e 


-N5e 


0 as 5 —oo. 


^ -2N Y. S 

'^USp{2N),qi^) ^ ^ 

DCA 
\D\<q 


PI I Z{D,D)Z{-D,-D)Y{D) 

^ ’ y y{-d)z^{d-,dy 


R 


Y X\HD{Wr 


A/D=WiVJ-VJWii r=l 
\Wr\<2 


Thus we can deduce 
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Theorem 7.2. Let <5 > 0 and indicate integration from 6 — ioo to 6 + ioo Take F(xi, ■ ■ ■ , Xn) 
a holomorphic function of n variables such that 

(7.30) F (xi, • • • ,Xn) = F (zbxi, • • • , ±Xn) . 


and assume F{xi, • • • , x„) = / ■ ■ ■ , decays rapidly in each variable in horizontal strips. 

With defined as in Equation \ 7.8[ let (^i, • • • ,(,n) = 0 if |.^i| + • • • + \fn\ > 2q — e for some 
e > 0. Then 


(7.31) 


L 


USp{2N) 

1 

(27ri)^ 




3ir-,]n=-oo 


(2Ar)l^l 


AULUA7={1,... ,n} 


/ JuSp(2N),qizK'J-ZL)F{iZi,--- ,iZn)dzi---dZr, 

■7(<5) J{-sF J{o) 


where ^{A),F[DiWr),Y, Z are defined as above. 


7.2.1. Special Case q = 1. Again, we particularly want to consider when q = 1. Then in the sum 
within Jq{A) in Theorem 17.21 the only valid term is that which arises from D = it). 


(7.32) 


(7.33) 


where 

(7.34) 


-2N Y. ^ 

JI{A)=Y^ (_l)I^L 

DCA 

\D\<1 


IZ{D, D)Z{-D, -D)Y{-D) 

Y{D)Z^{D-,DY 


R 


Y X{HD{Wr 


A/D=WiU-VJWRr=l 

\Wr\<2 

R 

Y 

A=WiU---UWii r=l 
\Wr\<2 


770 (Wr) = 


(2a) 


Wr = {a} 


(t) (« + /5) = {a,/3}. 


8. Example Calculations 

The purpose of this section is to explicitly calculate the 1 and 2-level density functions of zeroes of 
families of L-functions of elliptic curves to illuminate the n-level density calculations in Section 15.51 
We will compare our solution with a previous result when n = 1 |48j and to the limiting behaviour of 
random matrices when n = 2|16j to demonstrate that our results match up with current knowledge. 
The sums of the form 0 < d < X in this section are over fundamental discriminants. 
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8.1. Calculation of Sf. Let O- denote the path from —<5 — oo up to —S + oo and let denote 
the path from h — oo up to (5 + oo. Let / be a holomorphic function of one variable such that 


( 8 . 1 ) 

Then 

( 8 . 2 ) 


f(e) = f(-0). 


0<d<X j>0 

Xd{-M)LJE=i 


where is the height of the jth zero of the L-function Le{s, Xd)- Then following the calculations 
in Sections 15.51 


2(27ri)5f (/) ={/-/) ^ 

lc+ Jc-J 0 <(i<A 

Xd{-M)u}E=+i 
JM| 


+ 0!,Xd)fi-ia)da 


(8.3) 

(8.4) 

(8.5) 


The final step follows from Theorem 15.21 We can move the integration onto the imaginary axis as 
there are no poles when we only have a single variable. Hence 


E / / JEizK'd -ZL,ZM)f{iZl)dZi 

^^ InK Ir'LuM 

KULUM'''-'+ 

={ 1 } 

E / / Je{zk -ZL,ZM)f{iZl)dzi 


KULUM 
={ 1 } 

+ o(x^+^ 


( 8 . 6 ) 

(8.7) 


/ OD 

f{z) 

-OO 


d*E {^z, 0) + (—iz, 0) 


-J% {0,iz) 


dz + o[x^^^'^ , 


recalling that the definition of Jg {A,B) is 


' / — \ - E 2<5 

4{AB)= E 


0<d<X /3&B 

Xd{-M)‘^E=+^ 


DCA \ 


\ 2 tt I 


n 


r(l-^) Z^iD-,D-)Z^iD,D)Y^iD) 

r(i + 5)V' 


- V zliD-,DrY^iD-) 


R 


x{-ir\AE{D,D,D) 

AID= r=l 

WiU---UWe 


( 8 . 8 ) 
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where 


(8.9) 


( 8 . 10 ) 


( 8 . 11 ) 


Hd (Wr) 

Hd (Wr) 


Hd {Wr) + Ad^I (Wr) + Ad, 2 (Wr) + Ad, 3 (Wr) 

^SeD (1 + a — <5) — ^ (1 + a + 5)^ — ^ (1 + 2a) 

< ^ (1 + Ol + 0 , 2 ) 

1 

0 

\ 


d 


AD,^{Wr)= n ^^OgAE,i{A,B,D) 




B=A 


and Ae,Z(^,Y(^ etc defined as in Equations 15.111 - 15.141 


Wr = {a} 

Wr = {01,02} 
Wr = 0 
|1E^| > 3 


8.1.1. Checking Sf against Known Results. We want to compare this with the result from [48] 
where the authors use the ratios conjecture to calculate the one-level density for a family of elliptic 
curves. They show agreement numerically between their formula and the empirical one-level density 
of the zeros over a long range, in regimes where random matrix behaviour dominates and as the 
arithmetic terms take over in importance. 


Theorem 8.1 (Huynh, Keating and Snaith jUj). The 1-level density family of even quadratic 
twists of the L-function of an elliptic curve with prime conductor M is 


( 8 . 12 ) 


where 

(8.13) 


Siif) = 


1 


0 <d<A \ V / 

Xd{-M)u}E=l 


m E 


+ 2 


C L 

— — (1 -|- 2it) ^ isyinnf, 1 -|- 2it) {it, it) 

C Le 


( VM|d| ' 

V 2- , 

+ o(x^+^ 


-2it 


r(l — it)C{l + 2 it)LE{sym?, 1 — 2 it) 
r(l -|- it)LE{sym'^, 1) 


BE{—it, it) 


^ = E E Kvi) 

0<d<X j=—oo 
Xd(-M)ujE=+l 


for f, an even test function as described in (EH). ue is the sign from the functional equation of 
Le{s) as described in Section \4^ LE{sym?,s) is the symmetric square L-function associated to 
Le{s) and Be and B}, are arithmetic factors. The construction of the last three functions are 
described in [iS] (Note that Be, and B^ are called Ae and A}, in Huynh et al’s paper, but we have 
relabelled them to avoid confusion with our analytical terms.) 
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We notice that Si{f) differs from our definition of Sf{f) in (|8.2p . as Si{f) is a sum over all 
zeros of the L-functions while Sf{f) is only over the positive zeros. As we are restricting to even 

L-functions ( i.e Xd{—M)^E = +1)) we expect Si{f) = 25f'(/). 

This result matches perfectly with 5f(/). The term containing 2log ^ (1 + zt) + 

/ 

^ (1 — it) corresponds exactly to (0, it). If we consider the parts of Sf{f) corresponding to 
(it, 0) + {—it, 0) when D = (1) we get 


(8.14) 


(8.15) 


(8.16) 

(8.17) 


fiz) X] Aij (0,0,0) H(i,{{iz}) + Hinii-iz}) 


dz 


0<d<X 

Xd{-M)ujE=l 


j (1 + 2iz) - ^ (1 - 2i2;) 


0<d<X 
Xd{-M)u>E=i 


+ ^ ^D,j{{iz}) + AD,j{{-iz})dz 


i=i 


/ OO 

fiz) 2 

-OD rv ^ \r 


0<d<X 
Xd(-M)u)E=l 


(1 + 2i2;) + '^AD,j{{iz}) 


i=i 


dz. 


This matches up with the term 


(8.18) 


2 


C L 

(1 + 2iz) + —^ {sym^, 1 + 2i2;) + B\{iz, iz) 

C 


in (I8.12p . In our definition of the Ratios Conjecture in Section [5.11 we defined the analytical term 

to contain the term that Huynh et al. separated out from their analytic term 

Be- Checking the constructions of the Ratios Conjectures in their paper against ours, it is clear 

/ 

that ^DjUiz}) = ^ (sym^, 1 + 2i^) + ^^(iz, iz). 

All that remains is to check that the remaining terms in (I8.12p match up to the terms in 
J|;(iz,0) + J^(—iz,0) where D ^ tj). This is clear once we see from the construction of the 
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arithmetic terms that BE{—iz,iz) = AE{{iz} , {iz}, {iz}) . So 


(8.19) 


( 8 . 20 ) 


/ OD _ 

/(^) E 

-CO r, ^ 


^<d<X 
Xdi-M)LJE=i 




(VM\d\\ T{l + iz) o- /r • t r • i r - in 
I 27 r J r(l “ 2iz)Ae{{-iz}, {-iz}, {-iz}) 


dz 


/ OO 

/(‘) E 2 

■°° 0<d<X 

X(i( —= l 

LE{sym?, 1 — 2 i^) 

LE{sym?, 1 ) 


r(i - 2i^) 


\ 271 


BE{—iz, iz) 


r(i + iz) 


dz. 


Hence it is clear that our result in Theorem 15.111 for re = 1 agrees with the previous result in 
Theorem 18.11 


8 .2. Calculation of As an further example, we demonstrate the calculation of the 2-level 
density of zeroes of families of quadratic twist elliptic curve L-functions with even functional equa¬ 
tion: a family with orthogonal symmetry. Let C_ denote the path from —5 — oo up to —6 + oo and 
let C+ denote the path from 6 — oo up to <5 -|- oo. Let / be a holomorphic function of 2 variables, 
such that 

( 8 . 21 ) / ( 01 , 02 ) = /(± 01 , ± 02 ). 

Then 

(8.22) S^{f)= Y. E 

0<d<X ji,j2>0 

Xd{-M)ojE=+l 


where 'jj^d is the jth zero of the L-function LE{s,Xd)- We know from Cauchy’s Residue Theorem 
that 


(8.23) 

2\2Ei)^Si{f) 


'C+ 


IC- 


E 


0<d<X 
Xd( — 




/3,Xd)f{-ia, -iP)dadl3. 
















76 


A.M. MASON AND N.C. SNAITH 


We note that the functional equation of LE{s,Xd) can be arranged to give xj(l — 


^{s,Xd) - ^{s,Xd)- 

{2TrifS2if) 




(8.24) 


0<d<X 
Xd(-M)UJE=+1 

X / {—ia, —if3) dadj3 

-Jc L ^ ( 

JL.+ JL.- 0<d<X ' 

Xd{-M)uE=+l 

X / {—ia, —if3) dadfH 


Le \2 


V’ Le 


Oi,Xd 


Ltp ( 1 


Le \2 


+ /3, Xd 


Ic L ^ 


Le 


7- ./C+ o<d<X 

Xdi-M)ujE=+i 

X / (— za, —z/3) dadjd 


Le \2 


+ a,Xd 


■0 Le 


- l3,Xd 


+ 


LL s K-S 


0<d<X 
Xd{-M)ojE=+i 


2 - 


V' Le 


1 


- l3,Xd] f {-ia, -iP) dadjd. 


By the dehnition of Je{A,B) in (I5.92p . we can write this as 
2‘^i27TifSi^{f) 


'C+ 


{Je ({a, /3}, 0)) f{—ia, —i/3)dadj3 


(8.25) 


lc+ 


[ {Je {{/3,}, {-a}) - Je{{P, -a}, 0)) /(-*«, -i/3)dadf3 

Jc- 

f {Je ({a,}, {-/3}) - JE{{a, -/3}, 0)) f{-ia, -i(d)dadl3 


+ 


C- JC+ 

2 r 


<C- 


Je (0, {-a, -/3}) - Je ({-«}, {-/?}) 


-Je ({-/3}, {-a}) + Je ({-a, -/3}, 0) 


/ {—ia, —i/3) dadp. 


s,Xd) = 


Using Conjecture 15.21 we can replace the Je{A,B) with J*^{A,B) in all of the integrals above. 
In this section we regularly make use of the property that J)^{A,B) = J)^{A, B~): because of the 
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definition of ^(1/2 + P,Xd) in (15.661) the sign of the elements of B makes no difference. 

(8.26) ^ (-1)1^1 / [ J%{zk'J-zl,zm) f{-izi,-iz2)dzidz2. 

a:ulum={i,2} 

The next step is to move the contours onto the imaginary axis. We recall that /?}) has 

multiple poles as proved in Theorem 15.41 There is a pole at a = —j3 with residue ./|;({/3},0) + 
/3},0) — J|;(0, {/?}) and poles at a = 0 and /3 = 0 with residue 0. We then move these onto 
the real axis and combine similar terms, noting that f{x,y) = /(±x,±y) 

(8.27) 

22(27ri)25f(/) 

= ^^([ ) X] {-l)''^^''J*E{-iZK'diiZL-,-iZM)f{zi,Z2)dzidZ2 

\J-ooJ kuLUM 

={ 1 , 2 } 

/ + 00 

(J^({i/3}, 0) + J* ({-i/3}, 0) - 4 (0, {i/3})) f {13,13) dp + 0(XV2+^), 

-CO 

(8.28) 

= f[\ {-l)^^hE{-iZK^iZL,-iZM) f{zi,Z 2 )dzidZ 2 

VJ-oo/ kuLUM 

={ 1 , 2 } 

+ 22(27r)25f (/) + 0{X^/^+^). 

Given our definition of S^if) at the beginning of this section, it is clear 


(8.29) 

(8.30) 


Si{f)= Y 1 /(7ii,d,7j2,rf) 


0<d<X ji,j2>0 

Xd{-M}uiE=+l 


E Y.* In,d) + Sf if) 


0<d<X ji,j2>0 

Xd{ — ^)‘^E= + ^ 


where Yl* means terms with repeated variables are excluded. Define 


(8.31) 


SUf) = E E’ 


0<d<A ji,j2>0 

Xd(-M)ujE=+l 


This allows us to conclude that 
(8.32) 

1 


S*2{f) = 


2^{2'nf 


-oo/ KYdbYdU 

={ 1 , 2 } 


E (- 1 )'^'{-izK u izL, -izu) f{zi,Z2)dzidz2 + 0 (X^/^+") 


as shown in Theorem 15.111 
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8.2.1. Comparing S'® to known limiting behavior. Finally, we want to check that the limit of S^if) 
as X —oo is as expected i.e. matches that of the orthogonal random matrices from |16] . With 


(8.33) 


L = log 



consider 

(8.34) 


lim — 

X^oo X 


* 


E E'/ 


0<d<X ji,j2>0 

Xd(-M)ujE=+l 


’ I ’ 

TT ■^ / 


where is the imaginary part of the jth zero of LE{s,Xd) on the critical line and 


(8.35) 


From (j8.32jl we have 


V = E 1- 


0<d<X 

Xd{-^)‘^E=+i 


E E / 

0<d<X ji,j2>0 

Xd{-M)uJE=+l 

roo \ 2 


* f f 7il,d-h 7j2,d^ 


TT IT 


1 


IGvr^ 


([ ^ ^ {-l)^^^JEi-iZKCiZL,-iZM) 

\J — OO/ T^, I 7 - I ry-y 


(8.36) 


A'ULUM={1,2} 


X / ( Zi — ,Z2— ] dzidz2 + 0(X^/^+^) 


1 


TT TT 
2 


16L2 


(r\ V 

\J-QoJ ,r, ,jT_ri m \ L L L J 


A'ULUM={1,2} 

X / {ZI,Z2) dzidz2 + 0(X^/^+^). 

We evaluate the resulting sum making the following two large X approximations. Firstly that 


(8.37) 


7 E !(d) 


T)<d<X 


1 


E /(<«)■ 


T)<d<X 

Xd{-M)uE=+l 


Secondly, by the Euler-Maclaurin formula 
(8-38) 4 ^ ^ 


0<d<X 
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After the change of variables we will have a l/L^ outside the integral so we are only interested 
in the terms of order and higher, as the other terms will become negligible as X 


oo. 


n 


{ 0 <1 

f —nia 

—nib \ 

r’l 

1 i ’ 

L J 


(8.39) 


E 

0<d<X 

Xd(-A7)i^£=+1 


Ip f 1 nia \ /I nib 
\ 2 ~ ^ V2 ~ X’' 


‘Is —1 "p/3 




where ip {s,Xd) = Xd{—M)ujE i ^ is from the functional equation of LE{s,Xd) and 

so we recall 


(8.40) 


V’ A 

y (2 “I = 


, fVM\d\\ r', , r', 


We have 


(8.41) 


-Jt 


—nia —nib 


X* [ L ’ L 

If(—(") 4 »«> 4 o-«) 

4A. 


Next we note that for large X 


(8.42) 


^ 27ria\ L 

C IH-~- 

\ L J 2nia 


and 


(8.43) 


c \ L J 2nia 


and so consider terms of the form 
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(8.44) 


(8.45) 

(8.46) 



0<d<X 
Xd(-M)ojE = i 


(' ( 2ma\ 

-K'+— 



2'Kia\ , 

i + —)AE 


-Kia 


nib 

~T' 


Xd 



L 


2nia 


-dt 




nia 



L2 


~-(e 

ma 


— 27 ria 




where we recall that Ae{D) = Ae{D, D, D) and we have introduced AD{Wr) = A]^^i{Wr) + 
^D,2(14^r) + ^D,3(14^r') (See Theorem 15.21 for definitions). We also use the fact that Ae{D) = 1 if 
all the elements of D have value zero. 


So 

(8.47) 





, ^nsm(27ra) 

= -4L"^-^-4 

27ra 


Proceeding similarly with the other terms, and with the help of computer algebra, we can 
show that 


(8.48) 


lim —S^if) 
X^oo X* 


1 

4 



2 

/(6 *i, 02) det 7^50, ei;en d9id92 


where 

(8.49) 


Kso ,even {x,y) 


sm{n{y — x)) 
n{y - x) 


+ 


sin(7r(y + x)) 
n{y + x) 


This is exactly the form of the kernel for correlations of eigenvalues of random matrices from 
SO{2N) |16] and hence we have confirmed our result agrees with the conjectured limiting behavior 
of the family of zeroes [57]. 
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